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For  Peter  Jerome  Bellaire 

Stranger  and  Brother, 
Thanks  for  the  Clarinet. 


DIGEST 


For  an  aerosol  atmosphere,  disparities  in  gas 
temperature,  T  ,  and  particle  temperature ,  T  ,  may 

o  H 

arise  over  the  diurnal  cycle.  A  radiative- 
conductive-convective  heat  transfer  simulation  of 
the  dusty  atmosphere  of  Mars  was  developed  to 

quantify  the  temporal  and  spatial  character  of  the 

< 

temperature  difference  (T  -T  ) . 

H  O  . 

This  program  computed  thermal  energy  fluxes 
in  the  temporal  and  vertical  domains  in  a  51- 
layer,  vertically  inhomogeneous,  plane-parallel 
model  atmosphere,  with  a  30-layer  simulated  ground 
structure.  Observations  from  Mariner  and  Viking 
spacecraft  were  used  in  conjunction  with  other 
research  findings  to  provide  model  input  and  set 
the  boundary  and  initial  conditions  for  the 
simulation.  Temperature  fields  within  the 
atmosphere  were  inferred  from  radiative-conductive- 
convective  flux  fields  by  means  of  enthalpy  rate 
principles,  while  the  radiation  laws  of  Kirchhoff 
and  Planck  were  applied  to  the  flux  fields  to 
determine  aerosol  temperatures.  Several 
independent  methods  of  validation  were 
successfully  applied  to  the  model  output, 
including  comparisons  with  distinct  spacecraft 
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observations  and  the  separate  computation  of  dust 
heating  effects  outside  the  model.  <■ 

Just  after  local  noon,  dust  temperatures  at 
the  highest  altitudes  were  found  to  exceed  ambient 
temperatures  by  at  least  100  K,  while  in  the 
lowest  and  dustiest  atmospheric  layers, 
differences  were  typically  smaller  than  10  K.  The 
mean  temperature  of  these  lowest  and  thickest  dust 
layers  is  about  20  to  30  K  lower  than  surface 
temperature.  The  quantity  (Tp-T  )  decreased  with 
increasing  dust  load;  Tp  attained  a  minimum  near 
the  9  km  level,  with  the  altitude  of  this  minimum 
seemingly  independent  of  dust  load. 
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1.  INTRODUCTION 


1.1  Overview 

Certain  properties  of  gas-particle  mixes  in 
dust -laden  atmospheres  would  seem  to  be  uniformly 
applicable  over  a  wide  range  of  combinations.  The 
terrestrial  and  Martian  atmospheres  show  some 
similarities  in  the  character  of  their  respective 
gas-aerosol  interactions,  despite  the  obviously 
different  conditions  on  the  two  planets.  The 
primary  effects  of  atmospheric  aerosol  loadings  are 
radiative  in  nature,  and  this  thesis  examines  the 
manner  in  which  aerosols  can  modify  the  propagation 
of  radiation  in  model  atmospheres. 

Atmospheric  radiative  transfer  is  also  an 
essential  element  in  planetary  heat  distribution. 
The  processes  of  emission,  absorption,  and 
scattering  of  radiant  energy  are  different  for  the 
gaseous  (molecular)  and  nongaseous  (aerosol) 
components  of  an  atmosphere,  and  continuing 
research  concentrates  on  how  these  two  regimes 
interact.  Since  radiative  transfer  is  the  stuff  of 
which  climate  is  made,  this  approach  can  lead  to 
very  broad  conclusions  despite  its  limited  focus. 


2 


1.2  Statement  of  the  Problem 

Knowing  the  size  distribution  and  material 
composition  of  an  aerosol  leads  to  extinction  and 
scattering  functions  for  these  particles  via 
classical  Mie  theory,  whereas  knowing  the  molecular 
concentration  of  the  gases  involved  allows 
development  of  analogous  functions  for  these  gases 
using  Rayleigh  theory.  Quantification  of  local  and 
global  radiative  transfer  is  achieved  by  combining 
these  functions  with  multiple  scattering  effects, 
surface  albedo,  and  ground  emission 
characteristics . 

Radiative  fluxes  can  then  be  inferred  and 
enthalpy  (internal  energy)  balance  applied  to 
describe  the  spatial  and  temporal  distribution  of 
temperature  by  conductive,  convective,  radiative, 
and  phase  change  heating/cooling  processes.  Because 
the  various  atmospheric  components  have  different 
specific  heat  capacities,  deduced  heating  rates  for 
these  components  will  generally  not  be  the  same. 
Enthalpy  rate  disparities  arising  from  the 
differential  heating  of  the  aerosol  and  gas  act  tc 
promote  a  significant  terrperature  difference 
between  these  constituents. 


Research  has  shown  that  particle-gas 
temperature  differences  can  approach  values  over 
100  K  in  the  Earth's  mesosphere  (Fiocco  et  al., 
1975).  Mars,  having  an  atmosphere  that  is  very 
radiatively  active  (95^  C02),  may  exhibit 
temperature  differentials  quite  distinct  from  the 
terrestrial  case,  especially  since  its  aerosol  is 
persistent  and,  during  global  dust  storms,  attains 
much  higher  optical  depths  than  are  typical  on 
Earth  (Thorpe,  1981;  Zurek,  1981,  1982). 

1.3  Thesis  Objectives 

As  a  logical  extension  of  the  preceding 
discussion,  the  topical  objective  of  this  study  was 
to  determine  the  distinction  between  gas  and  dust 
temperature,  as  a  function  of  time  and  altitude, 
for  the  Martian  atmosphere  in  the  mature,  steady- 
state,  global  dust  storm  phase. 

The  methodological  objective  of  the  thesis  was 
to  construct  an  algorithm,  based  on  sound  physical 
principles,  to  compute  dust  temperatures  from 
output  fluxes  produced  by  the  MARSTEM5  model  of 
Failmann  (1976,  1977.  1983).  MARSTEMP  calculated 
solar  and  planetary  fluxes  as  a  function  of  time 
(over  the  diurnal  cycle)  and  vertical  coordinate 


(within  a  51-layer  plane-parallel  atmosphere  of  50 
km  thickness  and  a  50-layer  ground  structure  or  60 
cm  thickness) . 

Given  the  setting  of  a  tenuous  but  nearly 
chemically  homogeneous,  radiatively  active 
atmosphere,  justifiable  simplifications  to  the  dust 
temperature  algorithm  were  made  possible.  Direct 
observations  of  the  Martian  surface  and  atmosphere 
provided  by  the  Viking  and  Mariner  spacecraft 
permitted  reasonable  quantification  of  the  optical 
properties  of  the  dust  aerosol  using  Mie  theory. 

The  highly  ref ined  and  extensive  databases  of 
atmospheric  spectroscopic  line  parameters  compiled 
by  McClatchey  and  Rothman  provided  accurate  input 
for  the  simulation. 

An  additional  objective  of  this  investigation 
was  the  validation  of  model  findings  with  respect 
to  Viking  observations.  Model  output  agrees  rather 
well  with  conclusions  drawn  in  the  late  1970 's  by 
various  principal  investigators  immediately 
following  the  Viking  mission,  and  is  within  the 
numerical  range  of  actually  observed  or  reasonably 
extrapolated  values.  Validation  of  the  model  was 
also  achieved  by  several  independent  methods, 
including  the  comparison  of  dust  free  with  "dusty-' 
outputs,  which  were  mutually  consistent.  As  a 


2.  MODELING  CONCEPTUALIZATION  RASED  ON  REVIEW  Or 
THE  PERTINENT  LITERATURE 
2.1  Overview 

Modeling  radiative  transfer  for  a  vertically 
inhomogeneous  planetary  atmosphere  requires 
knowledge  of  the  radiative  output  of  the  parent 
star  (or  stars),  as  well  as  any  other  major-  energy 
sources  or  sinks  that  may  be  operating;  the 
chemical,  physical,  and  spectroscopic  properties  of 
the  planetary  atmosphere-ground;  subsystems;  and  the 
orbital  and  rotational  characteristics  of  the 
planet.  The  effects  of  aerosols,  variable 
atmospheric  and  surface  albedo,  olcanic 
outgassing,  clouds,  trace  gasas,  precipitation  and 
phase  changes,  loss  of  atmospheric  gas  to  space  or 
the  gain  of  extraplanetary  material  by  the 
atmosphere,  and  so  on,  need  to  be  accounted  for 
whenever  applicable.  A  comprehensive  model  would 
involve  a  multidisciplinary  effort  in  physics, 
chemistry,  astronomy,  astrophysics,  geology, 
geophysics,  meteorology,  etc.,  as  well  as  a  host  of 
related  specialties  arc  outfields.  Such  a  model 
remains  a  goal  as  yet  ur .attained. 

Ary  meaningful  analysis  requires  integration 
-T  ac .urate  data  with  sound  concepts  and  theory. 
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'ihe  physical  scale  of  the  phenomena  is  vast;  even 
the  most  recent  advances  in  technology  have  only 
made  it  possible  to  obtain  in  situ  observations  for 
our  nearest  planetary  neighbors,  and  the  Moon  alone 
has  received  human  explorers.  Acquisition  of  data 
for  planets  beyond  our  solar  system  is,  at  best, 
far  in  the  future. 

However,  given  the  chemical  and  physical 
nature  of  the  surface  (and  potential  surface- 
derived  aerosols)  of  any  planet,  the  scattering, 
absorption,  and  surface  albedo  characteristics  of 
such  material  can  be  ascertained  and  quantified  by 
current  theory  to  a  reasonable  degree  of  accuracy. 
Identification  of  the  constituent  gases  and  their 
abundance  by  volume  permits  the  determination  of  a 
mean  molecular  weight  for  the  atmosphere,  which  can 
be  used  with  the  equation  of  state  to  calculate  the 
thermodynamic  variables;  these  values  of  pressure, 
density,  and  temperature  indicate  directly  the 
local  effects  of  radiative  transfer,  and  also  give 
a  basis  for  constraining  the  model  to  a  particular 
regime  of  spectroscopic  interest  (strong-line,  weak- 
line,  or  intermediate-line?  Lorentz  profile, 

Doppler  profile,  or  Voigt  profile?)  for  atmospheric 
absorption  or  transmission.  Spatial  gradients  of 
the  pressure,  density,  and  temperature  fields  can 
also  be  derived,  shedding  light  on  the  fluid  flow 
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dynamics  that  create  global  circulation  and  climate 
patterns  that  ultimately  interact  with  radiative 
transfer  processes  in  a  nonlinear  feedback  mode. 

This  interaction  between  radiative  energy  and 
the  planetary  atmosphere-ground  subsystems  must  be 
a  paramount  consideration;  it  should  be  included  in 
any  comprehensive  model,  either  througn  an 
approximation  or  parameterization  of  the  feedback 
effects. 

2.2  The  Solar-Planetary  System 

The  sun  experiences  both  short-term  variations 
(as  evidenced  by  the  11\22  year  sunspot  cycle)  and 
long-term  variations  (as  evidenced  by  the  well- 
known  Maunder  Minimum,  from  1645  to  1715,  as  well 
as  other  similar  occurrences,  such  as  the  Sporer 
Minimum)  in  its  luminosity  and  energy  output.  These 
variations  are  not  well  understood,  yet  the  solar 
output  is  generally  considered  constant  at  3.827  x 

pk" 

10  W  (Pasachoff,  19(9). 

The  sun’s  electromagnetic  radiative  energy  is 
approximately  distributed  by  wavelength  X  as 
follows:  50$  in  the  infrared,  \ >0.7  /am;  -'’0$  m  the 
visible  spectrum,  0.4  *.m<  \<0.7  Am;  and  .0$  in  the 
ultraviolet,  A. <0.4  urn  (Liou,  1979).  In  these 
wavelength  regions,  the  sun  exhibits  differing 
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equivalent  blackbody  temperatures,  but  the 
effective  blackbody  temperature  is  taken  to  be 
approximately  5800  K.  Since  the  sun  is  considered 
here  mainly  as  an  energy  source,  its  size,  mass, 
and  other  physical  characteristics  will  not  be 
discussed  in  this  study. 

Mars  has  a  mean  distance  of  about  228  million 
km  from  the  sun,  with  a  revolution  period  (year)  of 
approximately  687  terrestrial  days.  Its  mean 
equatorial  radius  is  about  3397  km,  its  density 
about  3-9  x  10  kg/m  ,  and  its  mass  about  11%  of 
Earth's.  The  planet  has  an  oblateness  value  of 
0.009,  almost  three  times  the  value  for  Earth 
(0.0034).  One  Martian  day,  referred  to  as  a  "sol", 
consists  of  24  hours  37  minutes  or  88  775  seconds 
(terrestrial  units)  (Leovy  and  Zurek,  1979).  The 
planetary  surface  temperatures  range  from  about  140 
K  at  the  poles  to  probably  over  275  K  at  the 
equator  (Viking  measurements  at  the  surface  for  47° 
N  latitude  indicated  maximum  temperatures  over  240 
K,  and  for  22°  N,  maxima  near  260  K;  AGU,  1977). 

The  obliquity  (inclination  of  the  Martian 
equator  to  tne  Martian  orbital  plane)  is  23°59' , 
very  close  to  Earth's  value  of  23°27'.  The  orbit  of 
Mars  is  inclined  1°51'  to  the  ecliptic.  The 
eccentricity  of  the  orbit  (the  ratio  of  the 
distance  between  the  foci  of  the  elliptical  orbit 


and  the  major  axis  of  the  ellipse)  is  0.0934,  more 
than  five  times  the  value  for  Earth  (0.0167) 
(Pasachoff,  1979).  The  obliquity  varies  by  13° 
around  its  mean  value  over  time  scales  on  the  order 

5 

of  10  years,  while  the  eccentricity  varies  from 
near  zero  (almost  circular  orbit)  to  about  0.14 
(highly  elliptical  orbit)  over  time  scales  on  the 
order  of  10  ^  years.  Precession  of  the  Martian  axis 
occurs  at  a  frequency  nearly  in  resonance  with  the 
frequencies  of  several  of  the  planetary 
gravitational  perturbations  (mainly  Jovian) 
responsible  for  driving  the  obliquity  and 
eccentricity  variations,  thereby  creating  feedback 
effects  accounting  for  the  large  magnitude  of  these 
variations.  By  comparison,  the  amplitudes  of  the 
obliquity  and  eccentricity  variations  for  Earth's 
orbit  are  about  10  and  2.5  times  smaller, 
respectively,  than  lose  of  Mars  (Liou,  1979; 
Pollack  and  Toon,  1982). 

From  these  facts,  it  is  apparent  that  the 
Martian  orbit  is  far  more  irregular  than  Earth's, 
permitting  greater  extremes  in  insolation  changes 
and  seasonal  variations.  The  present  obliquity  and 
eccentricity  of  the  Martian  orbit  causes  the 
southern  hemispheric  summer/  northern  hemispheric 
winter  to  be  shorter  than  the  southern  hemispheric 
winter/  northern  hemispheric  summer.  Both  the 
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northern  hemispheric  winter  and  southern 
hemispheric  summer  are  warmer  than  the  same  seasons 
in  the  opposite  hemisphere,  since  the  former  occur 
at  perihelion  while  the  latter  occur  at  aphelion. 
Solar  radiation  influx  is  about  ^5%  greater  at 
perihelion  than  at  aphelion  (Pallmann,  1983). 

5 

The  range  of  orbital  deviations  over  10  to 

£ 

10°  years  allows  such  asymmetries  to  be  even  more 
pronounced  for  the  appropriate  combinations  of 
obliquity  and  eccentricity  extremes.  (In  this 
study,  tectonic  and  structural  changes  that  may 
occur  over  such  time  scales,  intrinsic  to  the 
planet,  and  their  effect  on  planetary  radiative 
transfer  characteristics,  will  not  be  considered.) 
Not  only  will  the  local  net  radiative  input 
oscillate  over  a  wide  range  under  such 
circumstances  (even  apart  from  any  intrinsic  solar 
variations,  which  we  have  already  neglected!),  but 
these  insolation  changes  in  turn  can  alter  the 
atmosphere-ground  subsystems,  which  then  may  affect 
the  net  insolation  in  a  complex  feedback  loop. 

As  a  consequence  of  such  irregular  orbital 
behavior  alone,  radiative  transfer  on  Mars  is  a 
highly  variable  and  nonlinear  process.  In  general, 
a  fully  time-dependent  radiative  transfer  model  for 
any  planetary  atmosphere  covering  characteristic 
time  scales  (even  if  inherent  planetary  structural 
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changes  are  ignored)  is  beyonu  the  scope  of  present 
knowledge.  Mars,  in  particular,  will  present  a 
challenge  for  generations  to  come.  On  human  time 
scales,  however,  these  changes  occur  so  slowly  that 
we  need  not  worry  that  our  models  are  restricted  to 
"present"  or  "constant"  climatological  conditions. 

2.3  The  Planetary  Atmosphere-Ground  Subsystems 

The  tenuous  atmosphere  of  Mars  is  largely  C02 

(95%  by  volume),  with  nitrogen  (2.7?)  and  argon 

(1.6$)  as  minor  constituents,  and  several  trace 

gases  in  negligible  amounts  (totalling  less  than  1$ 

by  volume)  (NASA/JPL,  198-4 ;  Hess,  1979).  This 

composition  leads  to  3  value  of  about  43 . 62  g/mol 

for  the  mean  molecular  weight  of  "Martian  air",  and 

a  value  for  the  specific  gas  constant  of  about  1.9 
2 

x  10  J/(kg  K).  A  surface  atmospheric  pressure  of 
7.7  hPa  (7.7  mb)  was  measured  by  the  Viking  1 
lander  at  22°  north  latitude,  but  at  both  Viking 
landing  sites,  the  pressure  was  observed  to  vary 
seasonally  by  as  much  as  29$,  due  to  the  growth  and 
shrinkage  of  the  polar  caps,  which  are  effective 
CQ2  reservoirs;  an  average  surface  pressure  around 
6.0  hPa  is  considered  typical  (NASA/JPL,  196--; 
Pollack  and  Toon,  198 2).  Both  CU2  and  H20  clouds 
have  been  observed  on  Mars  (Hunt  and  James,  1989  J . 


The  northern  ice  cap  (at  a  temperature  near  140  K) 
contains  a  permanent  water  ice  core  which  is  the 
main  source  of  atmospheric  water  vapor.  This  ice 
cap  allows  peak  values  of  50  to  100  precipitable 
M.m's  (column  abundance)  of  water  in  the  summer  near 
its  perimeter.  A  planetary  average  for  water  vapor 
abundance  is  taken  to  be  10  precipitable  yam's; 
during  simmer,  water  vapor  can  account  for  as  much 
as  two  to  three  percent  of  the  atmosphere  by  volume 
(Pollack  and  Toon,  1962;  NASA/JPL,  1984).  With  a 
high  degree  of  accuracy,  it  is  possible  to  restrict 
the  study  of  the  radiative  character  of  the  Martian 
atmosphere  to  the  study  of  C02,  the  most  abundant 
and  most  radiatively  active  constituent. 

Photographic  evidence  from  orbit  of  the 
erosional  action  of  liquid  water  supports  the 
conclusion  that  the  atmosphere  was  much  denser  on 
Mars  sometime  in  the  past  (Cordell,  1986). 
Calculations  indicate  that  the  atmosphere  has  lost 
about  10  times  as  much  nitrogen  and  about  20  times 
as  much  COP  as  is  present  today.  The  amount  of 
water  adsorbed  within  the  regolith  (weathered 
surface  material  or  "soli'')  and  existing  below  the 
surface  as  permafrost  is  unknown ,  but  considering 
this  and  the  permanent  north  polar  cap  of  water 
ice,  some  infer  that  Mars,  throughout  its  history, 
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may  have  pn  >du-  ••  *»i  enough  w  iter  l  eover-  the  entire 
surface  up  to  a  10  meter  aepth  (NAbA/JPL,  1984). 

The  surface  of  Mars  has  been  photogr'aphed  from 
orbit  b>  several  spacecraft,  and  in  situ  by  the 
Viking  landers,  whicn  also  subjected  soil  samples 
to  various  physical  and  chemical  tests.  The  reader 
is  referred  to  the  copious  literature  available  on 
the  Viking  missions  for  detailed  accounts  of  the 
results.  Grossly  oversimplified,  the  Martian 
surface  consists  of  a  highly  fragmented  layer,  up 
to  a  kilometer  thick,  of  boulders,  rocks,  dust,  and 
fine  particles,  all  probably  derived  from  volcanic 
basalts  (NASA/.JPL,  1984;  Pollack  and  Toon,  1982). 
Silicon  and  iron  are  the  most  abundant  elements; 
the  soil  is  about  45$  SiG2  and  about  19$  iron 
oxides.  Sulfur  is  about  100  times  more  abundant  in 
the  Martian  soil  than  in  Earth  soils,  but  potassium 
is  at  least  five  times  less  abundant  {NASA/JPL, 

19c-!  / .  The  finer  dust  particles  are  subject  to 
aeoiian  sail :  t  •  or.  md  lifting,  and  have  been 
observed  a*  n-  Igrr  ■  4J  kilometers  in  the 

atmosphere  during  '  hr  peri  odic  global  dust  storms; 
th*-?.-  parti  Te.  nr •  u ,  order  of  a  few  /urn’s  in 
radiu.-  ;  m.j>  :.»■  .->uf:  .  riently  rounded  by  aeoiian 
i  j.  over  1  .’.■-bit  y  "  j.v-uiogic";  time  scales  to 

t>-  c-  r i  .  ph-;  .  .  Pa  1  imann ,  1985),  although 

Vising  me.s  nr- merits  rj.jate  nonsphericity  is  more 
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the  norm  (Pollack  et  al.,  1979;  Chylek  and  Grams, 
1978).  The  atmosphere,  as  photographed  at  the 
surface  by  Viking,  has  a  reddish  tinge,  confirming 
that  the  haze  particles  responsible  are  indeed  Mie 
scatterers  (NSSDC,  personal  communication,  1986; 
NASA,  1978;  NASA/JPL,  1984).  Model  simulations  and 
direct  observations  of  Martian  dust  seem  to 
indicate  montmoril Ionite  219B  is  the  best  candidate 
to  approximate  its  mineralogical  and  optical 
properties  (Egan  et  al.,  1980;  Pollack  et  al., 

1979;  Martin  et  al. .  1979) . 

Periodic  global  dust  storms  have  been  observed 
on  Mars,  varying  in  severity  from  year  to  year. 
These  storms  originate  preferentially  in  the 
southern  hemisphere  during  the  summer  as  local 
storms  in  Noachis  and  Solis  Planum  (Peterfreund  and 
Kieffer,  1979;  PollacK  and  Toon,  1982;  Zurek, 

1982).  Gieraseh  and  Goody  .1973)  postulated  the 
"dusty  hurricane"  model  to  account  for  storm 
genesis  at  tnese  tropical  latitudes,  where  dust 
nesting  of  the  atmosphere  preyed  a  role  analogous 
to  latent  heating  in  terrestrial  hurricanes 
(Pallmann,  1976);  Viking  photography  has  indicated 
iocai  dust  devils  can  attain  heights  of  b  Km  or 
more  (Thomas  and  Gieraseh,  '.965).  in  a  few  week,.., 
these  local  storms  can  grow  to  global  proportions, 
aided  in  this  development  by  positive  feedback 


among  Una  dust  loading,  Uu-  .ism  n -i . it . -d  i non  1: .<  ■  1 
solar  heating  of  the  atmospher'e,  the  subsequent 
strengthening  of  the  storm  winds  and  mean 
meridional  circulation,  and  the  associated  thtrmal 
tides  (Pollack  and  Toon,  1962).  However,  the 
extreme  dust  loadings  that  are  ultimately  achieved 
(extinction  optical  depths  around  6  have  been 
observed;  Pollack  et  al.,  1979)  have  a  negative 
feedback  effect  on  the  global  storm  dynamics,  and 
the  storms  quickly  enter  a  plateau  phase  arid  slowly 
decay.  The  "e-folding"  time  for  dust  removal  (the 
time  required  for  a  b'j%  reduction)  from  a  20  km 
atmospheric  column  of  typical  dust  content  has  been 
calculated  to  be  about  56  days  (Haberle  et  al . , 
1982),  .implying  that  the  average  global  dust  storm 
has  a  lifetime  of  a  few  months  (on  the  order  of  a 
hundred  days  or  so).  This  agrees  with  observations, 
although  the  great  dust  storm  of  1971-1972  observed 
by  Mariner  9  was  notable  in  its  longevity.  However, 
at  altitudes  above  20  km,  dust  has  been  observed  to 
persist  throughout  the  year,  regardless  of  season. 

Intricate  layered  terrain  features  at  the 
polar  caps,  known  as  the  laminae,  so  far  have 
defied  adequate  explanation.  7rney  have  been 
theorized  to  be  dust  covered  mobile  water-ice 
sheets,  and  thus  may  be  intimately  related  to  the 
global  dust  storms  (Cordell,  198b).  Dust  from  such 


storms  is  preferentially  deposited  over  the  polar 
ice  caps  by  general  circulation  patterns  and 
through  the  scavenging  action  of  C02  condensation 
and  precipitation  from  the  atmosphere.  Dust 
deposition  has  significant  effects  on  the  ice  cap 
albedo,  which  again  brings  feedback  effects  into 
play  as  the  clearing  of  the  global  dust  storm 
accelerates  (Pollack  and  Toon,  1982). 

The  complex  relationship  between  insolation 
and  global  dust  storms  in  the  tenuous  yet 
radiatively  active  atmosphere  of  Mars  indicates  how 
radiative  forcing  and  planetary  subsystems  can 
interact  to  modify  each  other.  Meaningful  models  of 
radiative  transfer  should  account  for  such 
feedback,  even  if  only  in  parameterizations. 

2.4  Modeling  Radiative  Transfer 

Radiative  transfer  processes  involve  the 
emission,  absorption,  scattering,  and  transmission 
of  electromagnetic  wave  energy  in  a  given  medium. 

In  this  case,  the  modeling  of  radiative  transfer  on 
Mars  can  be  restricted  to  an  atmosphere  of  carbon 
dioxide  with  a  variable  optical  depth  caused  by 
seasonal  dust  loadings.  The  compilation  of  C02 
spectroscopic  absorptior.-line  parameters  by 
McCiatchey  et  ax.  in  1973  and  by  Rothman  et  al.  in 
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lytio.  Dei nnend jiari '  3  1969  i<  •rivation  of  i 
distributions  and  scattering  parameters  fur 
spherical  polydispersions,  the  radiative  transfer 
equation  evolved  earlier  by  Chandrasekhar,  van  de 
Hulst,  and  others,  and  the  well-known  algorithms 
for  single  and  multiple  ."lie  scattering,  are  the 
basic  tools  used  for  quantifying  radiative  transfer 
in  present  models  of  the  Martian  atmosphere. 
Continual  incorporation  of  new  facts  with  these 
tools  permits  attunement  of  the  models  to 
observations. 

Approximations  and  parameterizations  can  be 
used  in  the  models  to  reduce  computation  time  and 
effort,  and  may  actually  provide  insight  into  the 
physical  processes  involved.  Hie  determination  of 
directional  scattering  intensity  in  the  atmosphere 
via  phase  functions  that  either  require  truncations 
of  scattering  in  the  forward  direction  (Hansen, 

1969)  cr  approximations  using  the  delta-Eddington 
technique  (Joseph  et  al.,  1976;  Zurek,  1978)  has 
proved  advantageous  and  has  yielded  excellent 
results.  Indeed,  Hansen  (1969;  concluded  that 
optically  thick  atmospheres,  to  a  high  degree  of 
accuracy;  act  as  if  their  scattering  actually  is 
truncated  in  the  forward  direction,  due  to  the 
cumulative  effects  cf  multiple  scattering.  However, 
the  algorithms  for  multiple  scattering,  being 


rather  complex,  rieeo  improvement  with  respect  to 
describing  the  effects  of  the  variable  Martian  dust 
loadings.  Carbon  dioxide  absorption  lines  for  the 
tenuous  Martian  atmosphere  can  be  better  quantified 
using  the  Voigt  line-shape  profile  and  the  so- 
called  "intermediate-line"  approximation  (rather 
than  the  strong-  or  weak-line  exclusively) 

(Pallmann,  1976,  1983). 

Extension  of  well-known  terrestrial  layer 
models  (along  with  their  parameterizations 
associated  with  turbulence,  convection,  conduction, 
scale  height,  etc.)  to  the  Martian  setting  allows 
the  vertical  structure  of  the  atmosphere  to  be  more 
clearly  defined,  and  in  concert  with  thermodynamic 
equilibrium  considerations,  permits  the 
construction  of  vertical  temperature  soundings  for 
comparison  with  satellite  observations  (Pallmann, 
1976,  1977,  1983;  Egan  et  ai.,  1980).  This  is 
especially  important  when  attempting  to  describe 
dust  leading  effects  on  radiative  transfer.  The 
principles  of  terrestrial  aadley  cell  dynamics  and 
global  circulation  patterns  (Haberle  et  al.t  1982; 
Haber; e,  198o),  as  well  a..;  those  of  katabatic  winds 
union  tec  by  the  significant  meridional 
temperature  gradient;  w.  Fernandez ,  Saint  Louis 
University  MS  thesis,  ly'.'i)  ana  oust  hurricane 
formation  (with  dust  heating  playing  a  role 
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shapes  and  particle  corps.  .tiers,  Dot  at  present 


forms  the  core  of  undent,  -riding  for  particulate 
absorption  and  scattering.  A  better  mineral ogle 
match  to  Martian  dust  and  a  more  detailed 
spectroscopic  description  for  such  dust  would  help 
elucidate  the  complexities  of  the  radiative 
interaction  between  the  carbon  dioxide  atmosphere 
and  tno  <>.  usonai  dust  in  lection::. 

The  ’’scavenging”  of  aunt  from  the  Martian 
atmosphere  by  carbon  ilex  ice  precipitation  requires 
a  more  detailed  analysis,  with  emphasis  y.  the 
^hase  charge  durir g  deposition  aid  the  mechanics  of 
X:  ice  accretion  on  dust  nuclei .  The  albedo 
cnanges  of  the  polar  caps  a., res  by  preferential 
-ust  d' -position  t ;.<.<:•>  • . .•;<„•  r.ee  :  further  study, 

focusing,  atter.f  ;  or.  or  ‘  f-.-  a  .  o:  end  scat tr.cj.ng 
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of  C02  and  H20,  and  Line  nature  of  the  laminae. 
However,  for  the  purposes  of  this  study,  the  polar 
regions  and  their  associated  phase  changes  will  be 
left  to  future  study  and  the  model  applied  only  to 
mature  dust  storms  near  midlatitudes.  Changes  at 
midlatitudes  in  local  surface  albedo  due  to  the 
redistribution  of  dust  after  the  cessation  of 
global  storms  should  be  investigated,  since  these 
variations  could  be  significant  enough  to  affect 
local  thermal  equilibrium  and  alter  subsequent  dust 
storm  patterns  and/or  intensity  (Wells  et  al., 
1984);  this  issue,  too,  will  be  left  for  future 
consideration. 

2.5  Modeling  Convective  Heat  Transfer  on  Mars 

Pallmann  (1983)  describes  an  atmospheric  model 
for  Mars  in  which  an  interim  vertical  temperature 
profile  is  calculated,  then  suDjected  to  convective 
adjustment  by  determining  which,  if  any,  of  the 
model's  layers  in  this  temperature  profile  possess 
a  superadiabatic  lapse  rate.  Tee  enthalpy  of  these 
.ayers  is  redistributed  within  neighboring  layers 
above  and  below  until  all  such  layers  have  a 
strictly  adiabatic  lapse  rate.  The  layer  next  to 
the  ground,  however,  is  permitted  to  show  steep 
superadiabatic  lapse  rates,  in  accorcance  with  the 


accepted  concepts  of  Krai  chain's  turtuI-T.t  ttv 
convection  theory,  as  well  »s  KK.nirHjbukh.tv  length 
and  roughness  length  tn-'Or  y  v  Pa  I  imam,  I  Jo  -  ; . 

This  dual  approach  to  vertical  convective  heat 
exchange,  consisting  of  turbulent  therm..! 
convection  theory  and  the  simulated  ••••uistr.s  :t  it:, 
of  enthalpy,  soeriis  effective  in  modeling  thermal 
stratification  in  the  lower  Martian  atmosphere. 

2.6  Optical  Properties  of  Aerosols 

To  understano  the  effects  of  aerosols  on 
radiative  transfer  in  model  atmospheres,  two 
parameters  are  fundamental:  the  size  distribution 
and  the  complex  refractive  index  of  the  aerosol . 

Tne  size  distribution  (along  with  the  wavelength 
interval  under  investigation)  determines  the 
scattering  regime  to  be  modeled  (Rayleigh  or  Mie), 
while  the  refractive  index  describes  the  absorption 
and  scattering  cnaracter  of  the  aerosol  (Caale  and 
Grams,  1975).  In  addition,  the  crystalline 
structure  of  tne  aerosol’s  constituent  ’tatter  may 
be  of  consequence;  in  a  lithospheric  aerosol  of 
silicates,  for  example,  seme  of  the  incoming 
radiative  energy  can  be  distributed  tmorg  phonons, 
in  the  crystal  lattice  structure  (the 
"reststrahlen"  or  residual  ray  effect ) ,  a 


phenomenon  quite  distinct  from  true  absorption.  A 
silicate-type  aerosol  would  seem  to  have  a 
brightness  temperature  lower  than  expected  because 
the  energy  of  the  phonons  would  be  unavailable  for 
emission.  This  effect  is  apparent  for  terrestrial 
desert  sands,  and  must  be  taken  into  account  when 
attempting  to  derive  brightness  temperatures  of 
desert  surfaces.  Collisional  shock  effects 
experienced  by  the  parent  bodies  of  crystalline 
aerosols,  as  in  cratering  episodes,  could  affect 
the  lattice  structure  and  the  phonon  resonance 
energies  involved,  but  this  will  be  considered  of 
minor  consequence  for  our  purposes. 

For  the  case  where  the  wavelength  of  the 
radiation  under  study  is  sufficiently  small 
compared  to  the  radius  of  the  particle,  classical 
Mie  theory  describes  the  processes  of  scattering 
and  absorption,  as  formulated  by  van  de  Hulst, 
Sekera,  and  Deirmendjian,  Otherwise,  Rayleigh 
theory  applies.  The  sizes  of  the  lithospheric 
aerosols  considered  in  this  study  fall  into  the  Mie 
regime,  as  previously  discussed.  The  determination 
of  the  size  distribution  itself  depends  critically 
on  accurate  sampling  techniques  for  the  aerosol 
particles.  However,  the  remotely  measured 
characteristics  of  a  giver,  particle  population, 
such  as  optical  extinction,  albedo,  absorption, 


etc.,  can  be  compared  to  the  came  attributes  of 
well  studied  and  directly  measured  terrestrial 
aerosol  size  distributions,  fhis  permits  us  to 
apply  what  is  known,  about  terrestrial  aerosols  to 
the  cases  of  aerosols  in  planetary  atmospheres  in 
general.  In  the  case  of  Mars,  Viking  lander  imaging 
allowed  estimates  of  the  size  distribution  function 
to  be  made,  and  the  spectral  similarity  of 
montmoril Ionite  219B  to  the  dust  constrains  the 
dust's  complex  refractive  index  to  a  narrow  range. 

The  real  part  of  the  complex  refractive  index 
accounts  for  scattering  effects,  while  the 
imaginary  part  describes  absorption.  These  values 
are  typically  functions  of  the  wavelength  of 
radiation  under  consideration,  and  published  values 
for  a  given  substance  are  usually  statistically 
averaged  over  a  given  wavelength  or  frequency 
interval. 

Aerosol  particles  in  general  are  irregularly 
shaped,  and  this  departure  from  sphericity 
decreases  the  backscattered  radiation  (Liou  and 
Sasamori,  1975)-  Inhomogeneities  in  the  aerosol 
particles  would  also  be  expected  to  cause 
variations  in  scattering  and  absorption.  Models 
using  homogeneous  spherical  particles, 
incorporating  the  weli-<nown  Mie  scattering  theory, 
thus  give  upper  bound  vanes  for  backscattering  and 


absorption  (Liou  and  Sasamori,  1975;  Cadle  and 
Grams,  1975). 

2.7  Aerosols  on  Earth  ana  Mars 

For  the  lithospheric  aerosols  of  Earth  and 
Mars,  erosion  of  the  particles  by  saltation  and 
other  wind-generated  collisions  makes  the 
assumption  of  sphericity  and  the  application  of 
classical  Mie  theory  quite  reasonable,  although 
Pollack  et  al.  (1979)  supported  the  view  that 
smooth  platelike  particles  (typical  of  clays)  give 
a  better  fit  to  the  Viking  extinction  data.  Ine 
assumption  of  sphericity  may  still  be  only  a  small 
source  of  error  if  one  considers  the  random 
orientations  of  an  ensemble  of  such  particles. 

Cadle  and  Grams  (1975)  found  that  for  the 
terrestrial  aerosol,  the  number  of  particles  tended 
to  increase  rapidly  with  decreasing  particle  size 
for  particles  of  less  than  1  /an  radius.  This  could 
imply  that  on  Earth,  the  smallest  aerosol  particles 
are  the  most  important  radiatively  because  of  their 
sheer  numbers.  However,  increasing  particle  radius 
gives  a  rapidly  increasing  optical  cross-section, 
so  even  the  less  numerous  large  particles  could 
still  dominate  the  aerosol  optical  effects,  due  to 
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their  increased  scattering  and  absorption 


capabilities 


Cadle  and  Grains  (197b)  also  determined  that 


for  terrestrial  aerosol  size  distributions,  the 


greatest  contributions  to  scattering  and  absorption 


are  made  by  particles  with  radii  between  0,1  and 


2.0  /i.m.  Size  distributions  of  Martian  aerosols  seen. 


roughly  similar  to  Earth's,  so  these  conclusions 


could  apply  on  that  planet  also,  although  our 


model's  dust  temperature  algorithm  assumes  that 


within  the  highly  radiative  and  tenuous  atmosphere 


of  Mars,  particle  size  effects  are  rather 


negligible.  Pollack  et  al.  calculated  an 


"effective'1  particle  radius  for  Martian  aerosols  of 


2.5  Am  from  Viking  oats,  in  agreement  with  earlier 


conclusions  made  by  loon  et  al.  in  1977  from 


Mariner  9  IRIS  data  (Pol^acK  et  al.,  1979) 


However,  this  "effective"  radius  is  the  result  of 


an  analysis  of  a  broader  size  distribution. 


including  submicrometer  sizes,  with  a  modal  radius 


of  0.9  M m.  Many  other  authors  have  derived 


effective  .radii  closer  to  0.5  /on.  As  a  compromise, 


a  radius  of  1  a®  will  be  assumed  in  tnis  study. 


The  real  and  imaginary  parts  of  the  refractive 


ir.aex  determined  by  Lip.  et  al.  in  197  3  for 


terrestrial  soil  particles  (Liou  and  Sasamori, 


1975)  are  1.525  and  C.Co,  respectively,  with  a 


wavelength  around  0.5  Am.  on  In* ■  otn*  r  hand,  Ivlev 
and  Popova  11973)  calculated  "effective"  real  and 
imaginary  values  cf  l.t>9  and  0.02,  respectively,  in 
a  sophisticated  model  of  the  terrestrial  aerosol, 
over  a  wavelengtn  range  of  .a. 4  to  19  /nr,,  assuming 
dry  conditions  and  incorporating  a  wide  range  of 
molecular  species  in  the  aerosol .  A  perfect  match 
for  Martian  dust  has  yet  to  be  developed,  but  the 
results  obtained  by  using  montmcri 1  Ionite  219B  have 
proved  adequate.  The  results  of  investigations  by 
Singer  (1985)  indicate  that  a  certain  class  of 
Hawaiian  paiagonite  soils  ure  the  best  spectral 
match  to  Martian  sciis  yet  Known.  Paiagonite  is  an 
amorphous  hydrated  ferric-iron/siiica  gel  formed  by 
the  alteration  and  secondary  weathering  of  volcanic 
glasses  under  semi-aric  conditions,  and  may  be 
described  as  an  extremely  poorly  crystalline  clay. 
Singer  does  not  claim  paiagonite  as  the  universal 
Martian  soil,  however,  ar.d  stresses  that  more 
crystalline  ferric  oxides  may  exist  in  great 
abundance  locally.  Indeed,  Viking  Orbiter  data 
indicate  that  surface  regions  exist  better 
described  as  hematite-like  material.  It  is  hoped 
that  the  Mcirs  Geosciences/Cl imatology  Orbiter 
mission  scneauied  for  the  latter  part  of  t'ni3 
decade  will  shed  some  light  on  this  issue. 


The  existence  of  trn  terresti’ial  stratospheric 
aerosol  layer  has  been  known  since  its  discovery  by 
Junge  et  al.  in  1961  using  bai  ioon-borne  impactors 
vMcCormick  et  al.,  1982)  but  only  in  the  last 
decade  have  satellite-borne  experiments  been  used 
for  its  investigation.  Satellite  observations  have 
shown  that  meridional  circulation  patterns  carry 
aerosols  poleward  from  their  high  concentrations  at 
the  tropopause  in  the  equatorial  regions  and 
preferentially  deposit  than  near  the  poles. 

Ae:osols  can  act  as  nuclei  for  Polar  Stratospheric 
Clouds  iPSC's)  if  the  polar  stratospheric 
temperature  is  low  enough  ana  there  is  adsorbed  or 
otherwise  intermixed  water  in  the  aerosols 
(McCormick  et  al.,  1962 ) .  This  is  especially  true 
over  Antarctica,  with  its  strong  subsidence  cell 
and  extremely  low  temperatures .  This  polar 
deposition  of  aerosol  particulates  because  of 
Hadley  ceil  dynamics  cold  air  mass  subsidence 
is  not  unique  to  Barth .  Ma”s  aJso  exhibits  such 
scavenging  of  atrnosph*'  ri  dust  particles  am  the 
formation  of  high  aitUuue  clouds  over  its  pole-, . 
This  is  accentuate'",  by  arbor  dioxide  deposition  on 
the  dust  causes  by  too  extreme  cole  of  ’.he  Martian 
polar  atiaccphere ,  u‘ t imntely  .earing  to  remo/al  of 
the  dust  by  precipi*  /.’.or.,  to  w-,.r.  as  subsidence. 

It  is  possible  that  ar.y  planet  with  a  moderate 
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inclination  to  the  ecliptic,  ai.  adequate  rotation 


rate,  and  a  sufficient  meridional  temperature 


gradient  would  experience  a  similar  poleward  flow 


pattern  for  aerosol  circulation. 


Dust  storms  on  Mars  have  been  observed  from 


orbit  and  from  the  planet  surface  by  tne  Mariner 


and  Viking  missions.  Heating  of  th-.  atmosphere  in 


tne  vertical  by  dust  is  thought  to  be  an  important 


factor  in  the  evolution  and  maintenance  of  Martian 


storms.  Tne  presence  of  a  thick  dust  aerosol 


modifies  the  vertical  temperature  profile  during 


the  day,  increasing  static  stability  near  the 


surface  and  the  topmost  dust  layers  (Pallmann, 


1.976;,  as  well  as  enhancing  the  optical  extinction, 


low-level  atmospheric  heating,  and  surface  warming. 


Martin  ana  Kieffer  (1979:  noted  tnat  the  dust  load 


of  tne  1977  storms  creates,  optica,  depths  in  the 


"ange  i.5  tc  2.0,  compare a  to  a  median  value  of  0.3 


of  the  "clear'"  pre-storm  a  m^spnere .  From  the 


orbital  observations  of  the  liking  Infrared  Thermal 


Mapper  (IRTM)  a  the  it  /crn  v.^veiengtn ,  Martin  and 


Kieffer  deduced  that  po : ar  warmings  of  over  80  K  at 


about  25  km  cTtitud*-  ran  t>  <;  ;u.;eu  by  „he  presence 


of  the  c'ust  aerosol .  79.  s  through 


inversion  of  the  measur  >:• :  .  Am  rasia.  ce,  assuming 


an  ensemble  gray  body  err  sdvit.y  -.r  0.2  for  the 


dust  cloud  (this  is  nr.’ 


vl  9. a  part-id. 
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emissivity ! ; .  They  point..  <1  .;ut  that  this  p.  ,  cr 
warning  bears  striking  resemblances  to  t  sudden 
stratospheric  warmings  vSSW's)  on  Rarer.,  . r though 
the  causative  factors  are  undoubtedly  different. 

2.8  Radiative  Transfer  in  Aerosol  Atmospheres 

Radiative  transfer  in  an  atmosphere-ground 
system  involves  successive  absorption  and  re- 
radiation  at  various  levels.  Atmospheric  "counter 
radiation"  is  the  downward  re-emission  of  absorbed 
upweliing  planetary  emission,  which  in  turn  is 
dependent  on  the  amount  of  incoming  solar  radiation 
absorbed  at  the  planetary  surface.  This  primary 
absorption  is  affected  by  the  surface  reflectivity, 
or  albedo. 

Liou  and  Sasamori  (1975)  determined  that 
surface  and  aerosol  albedos  interact  to  yield  a 
"local  albedo",  ana  resultant  values  between  O.j 
and  0,4  constitute  a  transition  regime  for  the 
effective  atmospheric  heating  rates  caused  by- 
aerosol  haze.  For  local  albedo  values  below  0.3, 
hazy  and  cle  \r-  atmospheres  show  approximately  the 
same  degree  of  warming.  As  the  local  albedo 
increases,  however,  and  even  more  dramatically  as 
trie  value  exceeds  0.4,  absorption  across  the  solar 
spectrum  increases  in  a  hazy  atmosphere  as  opposed 
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to  the  clear  case.  Thus  the  hazy  atmosphere  gives 
rise  to  greater  total  atmospheric  warming,  although 
it  car.  promote  cooling  near  th-  surface.  Liou  and 
Sasamori ,  as  well  as  other  investigators,  have  also 
clearly  established  tout  absorption,  and  thus  total 
wanting,  increases  with  . erasing  solar  zenith 
angle.  It  should  be  jpL*n\m.t  that  local  warming 
effects  are  critically  aepenuent  or;  seasonal  and 
geographical  factors,  as  well  as  the  previously 
discussed  rr.iercphysical  and  optical  parameters. 

In  planetary  climate  studies,  too  absorption 
and  scattering  characteristics  of  the  atmosphere- 
ground  system  in  the  near  xR  wavelengths  (about  1 
to  10  mi)  are  of  primary  importance.  Generally,  the 
atmosphere  contributes  mo.t  strongly  to  the 
absorption  of  these  wavelengths  while  the 
lithosphere,  and  aerosols  derived  from  it,  dominate 
the  scattering  processes.  The  "equivalence  tneorem" 
of  van  de  Hu 1st  (1980)  slates  that  we  can  simplify 
the  problem  of  describing  a  gas-aerosol  mix  by 
ascribing  the  scattering  capabilities  cf  the  gas  to 
the  particulates  and  the  noorpti on  characteristics 
of  th*.-  aerosol  rt  the  go-.  y.,is  amo..r.tc  to  modeling 
trie  system  as  in  "equlvsi  .  v  mixture  of  an 
absorbing  nonocattering  gar  with  a  scattering 
r. on; bsorbir  g  aerosol 
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reducing  the  complexity  of  the  mathematical  models 
used. 

Models  of  atmospheres  with  high  optical 
depths,  caused  Dy  cloudiness  or  high  aerosol 
loadings,  should  properly  consider  multiple 
scattering  effects  tc  fully  describe  me  radiative 
transfer  processes,  and  this  can  complicate  the 
mathematics  considerably.  Parameterizations  such  as 
the  delta-Eddington  technique  (w'iscombe,  19Y7[a], 
1977Tb J ) ,  however,  seem  to  handle  this  problem 
quite  well,  without  any  significant  loss  of 
accuracy.  Also,  because  atmospheric  absorption 
exhibits  a  pressure  dependence,  the  Curtis-Goodson 
approximation  can  be  applied  to  simplify  the 
mathematical  treatment  of  vertical  inhomogeneity, 
allowing  models  to  take  the  vertical  variation  of 
pressure  into  account  by  substituting  multiple 
layers  of  homogeneous  stratification  for  the  actual 
pressure  profile  •.  Liou  and  SoosmorL ,  1975; 
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2.6  Modeling  the  Radi at:  on  and  neat  Budgets 

The  temperature.-  of  aerosol  part  spies  a 
planetary  atmosphere  d**por.  ::  ,n  tin-  energy 
exchanges  occur  i  rig.  Umov.  1  r .< •  5 » ;  :<*  ..s  •  r.e  mil  il:v- 
processes  of  emiss Lor.  and  ibsoi  otior.  for  t.i>. 


particles,  coilisionai  exchanges  with  the  ambient 
gas  (convectively  mixed  wherever  superadiabatic 
apse  rate?  occur} ,  and  the  exchanges  possible  due 
to  phase  transitions  of  the  ambient  gas  or  of  other 
volatiles  within  the  atmosphere  (Stephens,  1983; 
Mugnai  et  al.,  i >78}.  This  study  avoids  the  case  of 
p.:ise  change  effe.-ts  on  Kars .  For  the  terrestrial 
case,  r'iccco  et  al.  (1975,  1976)  calculated  that 
the  particle  temperature  can  exceed  the  ambient  gas 
temperature  by  as  much  as  LOO  K  in  the  mesosphere, 
er.d  this  temperature  differential  varies 
significantly  with  season ,  altitude,  and  time  of 
day.  The  rate  of  change  of  toe  particle  temperature 
defines  the  rate  of  enthalpy  change  for  the 
particle,  which  is  equal  to  the  net  divergence  of 
ail  me  thermal  fluxes  mentioned  above,  namely, 
phase  ;■ ha: go,  omission,  absorption,  and  coilisionai 
fluxes  I including  turbulent  convection) . 

lhc-  temperature  of  trio  ambient  gas  it 
Similar, j  dependent  or.  its  phase  changes, 
so  ..sort  .  >.-x  'hangers ,  a-id  one  emission  and 
us  o ”ptior.  ■ uuictet  of  its  constituent  molecules. 

;  ,f  tr.o  erx-  len.uty  differences  between 
* .  i:  y.  •  !.  ,,f  particles  -nd  ambient  gas,  the 

st-.*  os  of  local  trermobyoami  c  equilibrium  (LTE) 

•...  ;•■/<  i  by  each  separate  subsystem  arr>  not 
•  ;  K-'-nt, ,  :D;  <*  nigher  pressures  where 
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Fiocco  et  al .  (1975)  found  for  Earth  that  the 
temperature  differential  tended  to  increase  with 
altitude  and  decreasing  pressure,  then  diminished 
slightly  before  it  attained  a  constant  limiting 
value  where  energy  exchanges  between  gas  and  dust 
were  negligible,  a  characteristic  of  the  radiative 
equilibrium  of  interstellar  dust  (as  calculated  by 
van  de  Hulst  in  1946  and  1949). 

The  black  body  radiation  laws  are  implicit  to 
any  model  of  radiative  transfer  and  related  energy 
budgets.  To  apply  the  modification  of  a  gray  body 
assumption  to  these  basic  olack  body  laws  for  the 
aerosols  in  a  dust  cloud,  a  value  of  the  bulk  or 
ensemble  emissivity  is  needed,  a  much  different 
situation  from  determining  the  emissivity  of 
individual  particles.  This  can  be  a  significant 
stumbling  block,  especially  for  nonhomogeneous 
aerosols  with  poorly  Known  refractive  indices. 

Also,  as  was  noted  earlier,  phenomena  such  as  the 
residual  ray  effect  can  invalidate  black  or  gray 
body  assumptions  by  modifying  the  energy  transfers 
tnat  occur.  If  black  body  tneory  cannot  be  applied, 
modeling  radiative  transfer  becomes  virtually 
impossible. 

Obviously,  great  care  must  be  'oaken  when 
drawing  conclusions  from  satellite  observations  of 
radiative  processes.  The  emissions  of  dust  aloft. 


as  discussed  previously,  <-*011  id  yield  an  invert,  a, 
brightness  temperature  totally  different  from  that 
of  the  ambient  gas,  and  this  discrepancy  would 
become  even  more  severe  if  the  data  corresponded  to 
aerosols  at  increasingly  higher  altitudes.  If  the 
aerosol  was  highly  crystalline  in  nature,  the 
residual  ray  effect  could  put  even  the  inverted 
brightness  temperature  of  the  dust  itself  in  doubt. 
The  use  of  an  inaccurate  bulk  or  individual 
particle  emissivity,  as  mentioned  above,  would 
cause  similar  problems.  At  depth,  below  the  dust 
layer,  the  thicker  atmosphere  can  be  warmed  by  the 
dust  emissions,  and  the  fact  that  the  aerosol  is 
not  present  in  great  quantities  in  a  particular 
layer  should  not  lead  to  the  conclusion  that  the 
layer  is  free  of  the  aerosol's  influence.  This 
effect  is  mo  e  pronounced  for’  a  radiatively  active 
atmospheric  constituent  such  as  carbon  dioxide, 
because  of  its  enhanced  absorotion  capability  in 
the  infrared.  A  moderate  dust  layer  at  medium  or 
low  altitude  can  affect  the  lapse  rate  profile  of  a 
murh  thicker  atmospr.eric  layer  above  it,  and  modify 
the  static  stability  as  well.  The  implications  for 
convection  and  other  dynamic  el  ects  are  profound. 


\  ^  *-  --  *. 


3.  METHODOLOGY 
3.1  Overview 

For  an  aerosol  atmosphere,  disparities  in  gas 

temperature,  T  ,  and  particle  temperature,  T  ,  may 

materialize  over  the  diurnal  cycle.  A  radiative- 

conductive-convective  neat  transfer  simulation  of 

the  Martian  dusty  atmosphere  was  developed  to 

obtain  data  on  such  gas-aerosol  interactions.  This 

thesis  proposed  to  quantify  the  temporal  and 

spatial  character  of  the  temperature  differential 

(T  -T  ) ,  and  then  to  analyze  and  interpret  the 
P  S 

result  and  its  ramifications  on  the  atmospheric 
thermal  structure  of  Mars. 

The  use  of  well-known  radiative  transfer 
algorithms  within  the  simulation  established  some 
familiarity  with  meteorological  computer  modeling 
as  well  as  emphasized  the  physics  involved  in 
planetary  thermal  behavior.  These  algorithms  were 
chosen  on  the  basis  of  their  proven  efficacy,  sound 
theoretical  foundations,  and  basic  compatibility 
with  available  computing  equipment  in  terms  of 
financial  and  time  constraints. 

Such  algorithms,  modified  as  needed,  were 
utilized  to  construct  a  practical  computer  program 
to  produce  thermal  energy  fluxes  in  the  temporal 


and  vertical  domains  in  a  mod'll  of  th(>  Martian 
atmosphere.  Temperature  and  opacity  data  from 
Mariner  and  Viking  space  probes  were  used,  along 
with  orbital  and  physical  constants  given  by 
Souders  (1970)  and  Mechtly  (1973),  to  provide  model 
input  and  set  the  boundary  and  initial  conditions 
for  the  simulation.  Temperature  fields  within  the 
atmosphere  were  inferred  from  flux  fields  by  means 
of  enthalpy  rate  principles,  while  the  radiation 
laws  of  Kirchhoff  and  Planck  were  applied  to  the 
flux  fields  to  determine  aerosol  temperatures. 

3.2  Description  of  the  Model 

The  MARSTEMP  program  of  Pallmann  et  al., 
developed  from  1974  to  1976  under  NASA  Grant  NGR  26- 
006-042-2,  was  the  nucleus  of  the  current  model, 
MARSDUST.  The  reader  is  referred  to  the  works  of 
this  author  listed  in  the  Bibliography  if  an 
explicit  description  of  the  program  is  desired. 

Since  MARSTEMP  was  written  in  FORTRAN-H  and 
subsequent  modifications  were  not  entirely 
documented,  significant  effort  was  expended  in 
updating  and  debugging  the  program  to  run  under  the 
FORTRAN  77  compiler  now  in  use  on  the  department's 
MASSCOMP  computer.  Many  difficulties  also  had  to  be 
overcome  before  the  lengthy  but  necessary  HITRAN 


database  was  in  a  format  easily  usable.  Most  of 
this  work  preceded  the  development  of  the  algorithm 
to  compute  dust  temperatures. 

MARSTEMP  utilized  input  files  derived  from  a 
database  of  atmospheric  spectroscopic  parameters 
compiled  by  McClatchey  et  al.  at  the  Air  Force 
Cambridge  Research  Laboratory  (AFCRL)  in  1973  for 
use  in  the  their  LOWTRAN  code.  MARSDUST  similarly 
requires  data  for  C02  extracted  from  the  updated 
database  named  HITRAN,  developed  by  Rothman  et  al. 
(1983)  of  the  renamed  Air  Force  Geophysics 
Laboratory  (AFGL)  and  subject  to  periodic  revision. 
The  1986  HITRAN  version  used  in  this  effort  was 
kindly  provided  by  Dr.  Laurence  S.  Rothman  of  the 
Infrared  Physics  Branch  at  AFGL  to  the  Department 
of  Earth  and  Atmospheric  Sciences  at  Saint  Louis 
University. 

A  random  band  transmission  model,  based  on 
that  of  Goody  (1964),  was  developed  for  MARSTEM3 
utilizing  a  combined  Lorentz-Doppler  line  shape 
profile  in  a  treatment  of  intermediate  line 
intensities,  to  better  account  for  the  spectral 
gappiness  in  the  absorption  spectrum  of  the  tenuous 
atmosphere  of  Mars;  a  more  sophisticated  approach 
using  the  Voigt  profile  from  AFC-L's  FASCODE,  as 
described  by  Clough  et  al.  (1931),  proved  beyond 
the  3cope  of  this  study.  The  transmission  model 


used  parameters  extracted  from  the  McClatehey 
database,  namely,  line  intensity,  lower  state 
energy,  and  collision  half-width;  MARSDUST  also 
uses  the  same  quantities,  although  taken  from  the 
Rothman  HITRAN  compilation.  An  algorithm  for 
effective  absorption  by  an  assumed  height-depenaent 
polydisperse  silicate  aerosol  was  installed  in  the 
program  to  account  for  dust  effects  on  radiative 
transfer.  MARSTHAP  used  data  for  silicate 
polydispersions  calculated  by  Deirmendjian  (1969) 
in  its  treatment  of  dust  effects;  MARSDUST 
incorporates  data  derived  by  Viking  investigators 
wherever  possible  to  improve  on  MARSTEMP's 
parameterization  of  the  aerosol.  A  sophisticated 
model  for  multiple  scattering  was  beyond  the  scope 
of  available  computing  facilities;  a  decoupled  and 
modular  algorithm,  as  describee  by  Isaacs  et  al . 
(1986)  for  use  in  AFGL’s  FASCODE  and  LOWTRAN 
routines  (Clough  et  al.,  1961),  would  oe  a 
desirable  improvement  to  the  parameterized 
treatment  given  multiple  scattering  in  the  MARS7FMP 
and  MARSDUST  programs.  Such  an  algorithm  would  seem 
ideally  suited  for  accomodation  in  a  multitude  of 
radiative  transfer  problems ,  assuming  one  possessed 
equipment  with  the  computing  power  to  handle  the 
considerable  number  of  iterations  required  oy  the 


program. 


MARSTEMP  was  attuned  by  Pallmann  et  al.  to  the 
McClatchey  input  by  using  data  obtained  during  the 
period  November  1971  to  February  1972  by  the 
Infrared  Interferometric  Spectrometer  instrument 
(IRIS)  on  the  Mariner  9  spacecraft,  as  well  as  the 
IR  radiometric  and  S-band  occultation  experiments. 
Since  similar  IRIS  data  was  not  available  from  the 
Viking  mission,  attunement  of  MARSDUST  to  the 
revised  Rothman  database  was  achieved  by 
constraining  the  model  to  cyclical  diurnal  balance 
under  the  condition  of  no  dust  loading,  such  that 
temperatures  at  all  levels  at  24  hour  intervals 
were  approximately  unchanged.  The  new  attunement, 
although  opposite  in  sense,  is  comparable  in 
magnitude  to  the  prior  adjustment,  and  seams  quite 
reasonable  in  light  of  the  distinct  differences 
between  the  two  databases  used.  The  Rothman 
compilation  has  nearly  double  the  number  of 
absorption  lines  in  the  spectrum  as  the  McClatchey 
listing,  necessitating  the  opposing  trend  of  the 
current  attunement  process.  The  MARSTEMP 
adjustment,  in  essence,  was  a  calibration  implying 
a  wider  "effective"  line  half-width  and  an 
increased  "effective"  absorber  amount  due  to  the 
relatively  coarse  resolution  of  the  McClatchey 
database.  This  attunement  had  to  be  negated  and, 


indeed,  reversed,  in  accordance  with  the  improved 
spectral  resolution  of  the  new  HITRAN  database. 

3.3  Governing  Equations  and  Boundary  Conditions 

The  governing  equations  for  MARSTEMP  and 
MARSDUST  are  the  Equation  of  Transfer  (Liou,  1978; 
Siegel  and  Howell,  1971)  and  its  derivatives,  and 
the  Ehthalpy  Rate  Equation,  along  with  the  laws  of 
black  body  radiation,  spectroscopy,  "perfect" 
gases,  and  thermodynamics.  Within  either  model, 
enthalpic  heating  rates  are  equated  to  the  negative 
divergence  of  the  pertinent  net  fluxes  of  surface 
and  atmospheric  infrared  radiative  heating 
(amplified  by  a  non-radiative  term,  the  conductive 
heat  flux)  and  of  solar  radiative  heating.  These 
radiative  fluxes  are  obtaineu  tnrough  the  numerical 
integration  of  the  Equation  of  Transfer  over  the 
wavelength  range  0.2  yam  to  50  pm,  the  altitude 
range  0  to  50  km,  and  the  entire  solid  angle.  The 
Enthalpy  Rate  Equation  is  solved  numerically  at  52 
atmospheric  levels  between  0  and  50  km,  and  at  30 
subsurface  levels  between  G  and  60  cm.  Solutions  to 
these  equacions  are  ccflmon  in  the  scientific 
literature,  and  the  reader  is  referred  to  Pallmann 
(1976;  1977)  for  a  detailed  treatment  of  the 
mathematics.  It  is  important  to  note  the  nonlinear 


character  of  the  iritegru-di  f ferent La  1  equation 


describing  radiative  transfer.  No  analytical 
solutions  exist;  approximations  are  required  and 
feedback  effects  are  implicit. 

The  optical  characteristics  of  polydisperse 
aerosols  were  tabulated  by  Deirmendjian  in  1969, 
and  the  principal  investigators  of  the  Viking 
mission  have  used  his  modified  ganma  distribution 
to  describe  the  size  distribution  of  the  Martian 
dust  particles.  Analysis  of  photographic  data  from 
the  Viking  Lanaer  imaging  cameras  established 
reasonably  accurate  values  for  the  necessary  input 
parameters  to  the  gamma  distribution  function, 
which  allowed  Mie  theory  to  be  applied  toward  the 
compulation  of  values  for  the  optical  constants  of 
the  dust  ^PoliacK  et  al.,  1979).  Other  researchers 
used  radiance  data  obtained  from  the  Infrared 
Thermal  Mapper  (IRTM)  on  the  Viking  Orbiter  as 
boundary  conditions  to  define  the  optical 
properties  of  the  planet  as  a  whole  under  varying 
conditions.  Kieffer  et  al.  (1976),  using  the  20  /wn 
band  of  the  Viking  Orbiter  IRTM,  unplied  a  mean 
value  of  0.285  for  the  surface  albedo  near  the 
Viking  1  landing  site  (Pollack  et  al . ,  1979)  and 
implied  values  between  0.90  and  1.0  for  surface 
infrared  emissivity.  Pollack  et  al. ,  using  the 
Viking  Lander  imaging,  derived  explicitly  a  value 
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of  0.86  for  the  single  scattering  albedo  of  the 
dust  averaged  over  the  solar  spectrum  centered  at 
0.67  Am,  which  amounts  to  setting  the  average  dust 
absorptivity  at  0.14  in  this  spectral  region.  Hunt 
(1979)  derived  ensemble  infrared  emissivities 
between  0.4  and  0.7  for  dust  clouds  under  variable 
dust  loadings  and  zenith  angles,  implying  high 
infrared  emissivity  values,  perhaps  as  high  as  0.8 
to  0.9,  for  individual  dust  particles.  Thus,  by 
analogy  to  terrestrial  silicate  sands,  infrared 
emissivity  values  between  0.8  and  0.9  may  be 
considered  typical  of  individual  Martian  dust 
particles.  These  parameters  were  required  by 
MARSDUST  as  a  prerequisite  to  the  establishment  of 
an  algorithm  to  calculate  dust  temperature. 
Conversion  of  the  dust  extinction  coefficients  used 
in  the  MARSDUST  simulations  to  actual  optical  depth 
values  is  dependent  on  accurate  figures  for  the 
integrated  column  density  of  the  dust  as  well  as 
the  mean  radius  of  the  particles  (Pallmann,  1976). 
In  the  abscence  of  agreement  on  these  quantities, 
precise  values  of  optical  depth  will  not  be  quoted 
here;  observations  indicate  that  the  dust 
extinction  coefficients  used  in  this  study  are  well 
within  the  range  of  optical  depths  and  dust 
distributions  possible  on  Mars.  In  any  case,  trends 
in  atmospheric  thermal  behavior  can  indeed  be 


ascertained  from  the  data,  and  further  analysis  or 
refinement  by  others  is  welcome  and  encouraged. 

3.4  Initial  Atmospheric  Conditions  of  the  Model 

An  initial  condition  vertical  temperature 
profile  for  MARSDUST  was  deduced  from  the  Viking 
data  of  Seiff  and  Kirk  (1977)  and  the  initial 
profile  used  by  Pallmann  et  ai.  in  MARSTEMP. 
Temperatures  above  the  35  km  level  in  the  Seiff  and 
Kirk  profile  disagree  with  the  Pallmann  profile 
based  on  Mariner  9  data,  and  in  one  review  by 
Pollack  et  al.  in  1979,  no  values  above  this  level 
were  quoted  at  alll  MARSDUST  output  was  very 
unstable  at  the  upper  altitudes  when  the  Seiff  and 
Kirk  profile  was  used,  which  had  the  temperature  at 
the  50  km  level  at  an  usually  cold  and  physically 
suspect  148  K.  This  is  very  near  the  frost  point  of 
002,  where  the  transition  from  the  gaseous  to  the 
solid  phase  may  occur;  no  such  behavior  was  ever 
observed  by  Mariner  or  Viking  spacecraft.  In  the 
upper  layers  of  the  Seiff  and  Kirk  profile,  diurnal 
balance  proved  unattainable  and  intense  warming  of 
these  highest  layers  could  not  be  eliminated  in  the 
model  output.  The  Seiff  and  Kirk  profile  is  perhaps 
more  properly  a  trajectory  profile  rather  than  a 
vertical  one;  its  uppermost  temperatures  seem  to  be 
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more  characteristic  of  the  polar  region  or  at  least 
the  higher  latitudes,  which  may  indeed  have  been 
the  region  of  entry  for  the  deorbiting  lander,  but 
which  is  very  far  removed  from  the  final  landing 
site  located  at  subtropical  latitudes.  Because  of 
such  practical  concerns,  the  profile  used  in 
MARSDUST  was  a  hybrid  of  the  Seiff  and  Kirk  data 
and  the  Mariner  9  data.  The  Mariner  9  profile 
implies  upper  layer  temperatures  warmer  than  Seiff 
and  Kirk,  including  a  50  km  value  near  160  K;  this 
profile  was  adopted  by  Pallmann  for  MARSTEM3. 
MARSDUST  assumes  the  lowest  level  temperatures 
quoted  by  Seiff  and  Kirk  are  accurate  but  adopts 
the  original  values  of  Pallmann  (1976)  for  the 
uppermost  levels. 

After  initialization,  the  model  requires  that 
several  hundred  timesteps  be  run  before  numerical 
stability  is  attained  by  the  output.  Planetary 
fluxes  in  the  model  output  increased  slightly  but 
steadily  ove^  the  first  few  days  of  simulation,  as 
did  derived  temperatures,  but  then  temperatures 
were  observed  to  oscillate  mildly  about  median 
values  after  approximately  three  simulated  days.  It 
appears  the  model  is  more  sensitive  to  variations 
in  solar  radiative  input  than  planetary  output; 
derived  temperatures  after  three  days  varied 
typically  by  less  than  five  percent  from  their 


initial  values,  while  planetary  fluxes  increased  by 
about  ten  percent  over  the  same  time.  These 
disparities  were  lessened  by  having  no  dust  load 
and  aggravated  by  increased  dust  loading,  and  are 
apparently  a  measure  of  the  success  of  the  applied 
attunement  process  described  previously,  as  well  as 
the  affects  of  dust  loads.  In  the  discussion  of 
results  offered  here,  values  quoted  are  for  the 
fourth  day  of  simulation  or  later,  when  model 
output  is  assumed  to  be  more  or  less 
computationally  stable. 

3.5  The  Dust  Temperature  Algorithm 

MARSDUST  includes  an  algorithm  to  calculate 
dust  particle  temperature,  derived  from  principles 
set  forth  by  J.R.  Howell  and  R.  Siegel  (1969).  An 
explicit  formulation  of  the  dust  temperature  module 
is  given  in  Figure  1.  The  model  consists  of 
idealized  concentric  spheres,  where  the  net  influx 
of  energy,  Q^,  tc  the  particle  is  the  sum  of  the 
contributions  of  longwave  and  shortwave  absorption, 
integrated  over  all  incident  zenith  angles  with 
azimuthal  independence  assumed.  Because  of  the 
extreme  minuteness  of  the  particles,  instantaneous 
thermal  equilibrium  can  be  considered  a  good 
approximation,  such  that  Q  =  -Q^  .  Absorptivity 


of  0.14  for  solar  (shortwave)  flax  is  assumed  in 
accordance  with  Pollack  et  al.  (1979).  By  invoking 
Kirchhoff's  principle  that  a  good  emitter  is  a  good 
absorber,  the  absorptivity  is  set  equal  to  the 
emissivity  for  the  planetary  (longwave)  spectrum. 
For  all  data  presented  here,  an  emissivity  of  0.80 
is  assumed  in  analogy  to  terrestrial  desert  sands, 
and  by  considering  the  suppressing  effect  of  the 
silicate  lattice  structure  on  emission 
(reststrahlen  effect),  which  would  tend  to  lower 
the  emissivity  value.  As  a  sink  mechanism,  C02 
conductivity  of  heat  away  from  the  particle  was 
assumed  to  follow  a  Gaussian  distribution  function 
(height-square  logarithmic  profile),  as  described 
by  Pallmann  (1983).  Such  a  vertical  dependence  was 
justified  by  considering  the  nature  of  the  eddy 
exchange  processes  and  turbulent  mixing  responsible 
for  carrying  the  dust  aloft.  It  was  assumed  that 
the  same  physical  processes  that  allow  dust 
suspension  permit  the  effective  conduction  of  heat 
away  from  the  particle  by  the  atmosphere,  such  that 
smaller  dust  loadings  at  any  level  corresponded  to 
less  effective  C02  conductivity  at  that  level,  and 
vice  versa.  As  a  consequence,  dust  and  atmospheric 
temperatures  at  the  lowest  levels  are  always  equal, 
with  any  dust  heat  excess  quickly  being  lost  to  the 
atmosphere. 
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In  the  model ,  a  dust  particle  with  effective 

spnerical  radius  r  lies  within  a  spherical  gaseous 

shell  of  radius  R  .  The  value  of  R  is  physically 
g  g 

limited  to  be  never  less  than  r^  +  d,  where  d  is  on 
the  oraer  of  the  mean  free  path  of  a  gas  molecule. 
For  the  tenuous  atmosphere  of  Mars,  this  strongly 
suppresses  any  direct  thermal  interaction  between 

the  gas  and  the  particle,  since  R  is  nominally  an 

g 

order  of  magnitude  larger  than  typical  r  's,  even 
at  the  lowest  altitudes.  The  algorithm  does  permit 
dust-gas  feedback,  such  as  that  described  by  Fiocco 
et  al .  (1975),  but  only  when  R^  and  Tp  are  of  the 
same  order  of  magnitude,  or  when  the  gas  has  a  very 
low  emissivity  in  the  thermal  infrared;  for  the 
Mars  case,  the  surface  area  of  the  gaseous  shell  is 
always  much  larger  than  tne  surface  area  of  the 
particle,  and  the  term  involving  gas  emissivity  in 
the  formula  given  in  Figure  1  is  negligible.  It  was 
therefore  deemed  appropriate  to  neglect  this  effect 
in  the  MARSDUST  model,  considering  the  assignment 
of  a  typical  1  /im  radius,  implicit  within  MARSTEMP, 
tc  the  individual  oust  particle.  (A  flowchart  of 
the  MARSDUST  program  is  presented  in  Figure  2  to 
assist  the  reader  in  visualizing  the  various 
modules  and  algorithms). 


It  is  important  to  note  that  this  assumption 
does  not  imply  any  lack  of  synergy  between  gas  and 
dust  in  the  model  when  one  considers  overall 
atmospheric  warming.  Indeed,  this  thesis  supports 
the  contention  that  dust  contributes  to  warming  of 
the  atmosphere;  but  the  effect  is  primarily  the 
result  of  increased  atmospheric  absorption  due  to 
multiple  scattering  of  radiation  by  the  dust,  not 
by  the  heated  particles  directly  wanning  the  gas  by 
radiation  or  conduction. 

Radiative  heating  of  the  gas  by  the  dust  can 
indeed  occur,  and  may  become  significant  for  the 
higher  dust  concentrations,  but  this  feedback 
effect  is  neglected  in  the  current  model.  MARSDUST 
output  indicates  that  particle  fluxes  are  much 
smaller  than  solar  fluxes  and  usually  smaller  than 
planetary  fluxes,  which  are  the  two  most  important 
thermal  forcing  mechanisms.  Flux  errors  introduced 
by  this  neglected  effect  are  probably  on  the  order 
of  10$,  roughly  the  ratio  between  solar  fluxes  and 
atmospheric  layer  fluxes;  errors  in  temperatures 
derived  from  the  Stefan-Boltzmann  Law  would  be  even 
less  because  of  their  fourth-power  relationship  to 
fluxes.  It  is  hoped,  however,  that  future  work  will 
improve  on  the  current  algorithm  in  treating  these 


nonlinear  effects. 


4.  DISCUSSION  OF  RESULTS 


4.1  Overview 

The  reader  is  referred  to  the  Appendix  for  a 
comprehensive  graphic  presentation  of  the  results 
of  the  MARSDUST  simulations.  Please  note  that  LMST 
is  an  abbreviation  for  Local  Martian  Solar  Time, 
while  DUSTAMP  is  shorthand  for  Dust  Amplitude,  a 
measure  of  the  atmospheric  dust  loading.  LMST  0.000 
indicates  the  beginning  of  a  diurnal  cycle,  while 
LMST  24.000  signifies  the  end  of  the  same  24  hour 
period.  All  plotted  data  for  dust  correspond  to  a 
dust  emissivity  of  0.80;  output  utilizing  this 
emissivity  differs  by  less  than  3 %  from  output 
using  an  emissivity  of  0.90.  A  sensitivity  study 
showed  that  derived  dust  temperatures  were  only  a 
weak  function  of  emissivity  for  values  above  0.70 
and  strong  dependence  was  seen  only  for  highly 
unrealistic  emissivities  below  0.30.  The  values  of 
emissivity  referred  to  here  are  for  the  individual 
dust  particle  and  are  quite  distinct  from  the  bulk 
emissivity  values  quoted  in  the  literature  for  an 
ensemble  of  particles  in  an  optically  thick  cloud. 

For  a  simulation  without  dust  loading  at  all 
timesteps,  atmospheric  temperature  soundings  above 
9  km  showed  the  smooth  exponential  decline  with 
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height  characteristic  of  radiative  cooling.  The 
lowest  9  km  were  dominated  by  conductive  and 
convective  processes,  as  exhibited  by  persistent 
nearly  adiabatic  lapse  rates  (+5 .2  K/km)  and  strong 
nightime  inversions  in  this  region  (Figures  3  and 
4).  At  night  in  the  lowest  2  km,  a  temperature 
inversion  is  a  constant  feature  because  of  the 
rapid  cooling  of  the  surface  after  sunset  with 
respect  to  the  overlying  atmosphere.  This  inversion 
intensifies  to  lapse  rates  over  -17  K/km  just 
before  dawn,  and  does  not  completely  disappear 
until  well  into  mid-morning,  near  10.000  LMST. 
Without  dust  loading,  maximum  surface  temperature 
is  around  250  K  just  after  local  noon,  while  the 
minimum  just  before  dawn  approaches  185  K.  The 
inversion  in  the  lowest  layers  sets  in  very  quickly 
after  sunset. 

Dusty  atmosphere  simulations  were  run  with 

three  extinction  coefficients  (representing  three 

,  -7 

distinct  dust  loadings;:  1.2  x  id  per  cm,  3.0  x 

-7  -7 

10  per  cm,  and  6.0  x  10  per  cm  (Figure  5). 

These  correspond  to  the  moderate  to  heavy  dust 

storms  observed  on  Mars  by  the  Mariner  and  Viking 

spacecraft.  The  dust  load  in  eacn  case  was  set  to 

zero  at  an  altitude  of  42  km,  in  accordance  with 

the  observation  that  Martian  dust  loadings  rarely 

exceed  this  level.  This  assumption  creates  an 
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abrupt  discontinuity  in  the  atmospheric  temperature 
profile  at  the  42  km  altitude.  The  reader  is 
cautioned  thar.  the  temperature  values  in  the 
uppermost  layers  of  the  dusty  atmosphere  soundings 
are  only  as  realistic  as  the  previously  stated 
assumption.  It  does  seem  reasonable  that  such  an 
abrupt  change  should  occur  at  least  somewhere  in 
the  upper  atmosphere,  since  purely  physical 
considerations  imply  a  shallow  transition  region 
above  which  aerosols  of  a  given  effective  radius 
cannot  be  kept  aloft. 

The  shapes  of  the  temperature  soundings  for 

the  various  dust  loads  are  quite  similar.  However, 

with  increasing  dust  load,  layers  above  25  km 

appear  to  be  cooled,  while  the  atmospheric  layers 

below  25  km,  and  the  surface  itself,  show  a 

significant  warming  trend  (Figures  6-9).  A  model 

sensitivity  study  indicated  that  surface  cooling  is 

achieved  only  after  a  drastic  increase  in  dust 

concentration,  corresponding  to  a  DUSTAMP  of  1.2  x 
-4  -1 

10  cm  .  This  is  a  thousand-fold  increase  in 
particle  number  density  from  our  lowest  previously 
stated  DUSTAMP  value,  and  represents  a  dusr  load  of 
about  5000  particles  per  cubic  centimeter,  assuming 
an  effective  radius  of  one  micrometer.  Predawn 
minimum  surface  temperature  for  the  lightest  dust 
load  is  near  190  K,  while  for  the  heaviest  dusr 
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load,  Lt  Ls  about  200  K.  Predawn  temperature  at  the.- 
40  km  level  for  the  lightest  dust  load  is  about  17b 
K,  while  for  the  heaviest  dust  load  it  is  closer  to 
172  K,  At  LMST  13.900,  the  maximum  surface 
temperature  for  the  lightest  dust  load  is  253  K 
while  it  exceeds  260  K  for  the  heaviest  dust  load. 
The  effect  of  the  thicker  dust  loading  seems  to  be 
mainly  felt  near  the  surface,  where  the  dust 
contributes  to  strong  warming  but  also  promotes 
cooling  aloft  by  inhibiting  escape  of  planetary 
emissions  to  higher  altitudes.  The  increased 
multiple  scattering  induced  by  the  higher  dust  load 
enhances  absorption  by  carbon  dioxide  in  the  lowest 
levels  at  the  expense  of  the  upper  atmosphere.  By 
adding  to  the  atmospheric  counter  radiation  in  the 
lower  layers,  the  thicker  dust  loads  can  help  warm 
the  surface  throughout  the  diurnal  cycle. 

Even  a  cursory  inspection  of  the  dust 
temperature  profiles  (Figures  10  and  11)  shows  the 
swift  response  of  dust  temperature  to  the  influx  of 
solar  radiation,  regardless  of  dust  loading. 
Similarly,  the  surface  temperature  experiences 
drastic  shifts  between  nighttime  and  daytime 
conditions..  The  response  of  the  atmosphere  as  a 
whole  to  solar  input,  whether  or  not  dust  is 
present,  is  far  more  sluggish.  This  difference  in 
response  indicates  the  disparity  between  the 
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specific  heat  capacities  and  thermal  inertias  of 
the  solid  and  gaseous  components  of  the  planetary 
ground-atmosphere  system. 

With  increasing  dust  load,  dust  temperatures 
attain  lower  maximum  and  higher  minimum  values,  and 
th<_-  magnitude  of  the  dust-atmosphere  temperature 
differential  decreases.  The  higher  dust  content 
also  aids  in  burning  off  the  morning  inversion  more 
quickly.  Higher  dust  amplitudes  tend  to  lower  the 
altitude  of  minimum  dust  temperature  at  night,  and 
raise  it  during  the  day.  This  implies  that  the 
upper  levels  "feed"  off  the  heat  content  of  the 
very  dusty  lower  levels  throughout  the  diurnal 
cycle;  this  is  more  evident  at  night,  of  course, 
when  the  heat  source  is  not  replenished,  while  in 
the  daytime,  the  heat  content  of  the  lowest  layers 
steadily  increases  with  continued  insolation.  While 
dust-free  and  light  dust  load  conditions  give  very 
stable,  almost  isoth  ,  '.-ipse  rates  at  the 
topmost  levels,  increa  v-ju  dust  oading  yields 
increasingly  stable  lapse  rates  at  nearly  all  other 
levels  higher  than  9  km  or-  so,  throughout  the 
diurnal  cycle.  The  lapse  rates  in  the  lowest  9  km 
do  not  seem  to  be  a  strong  function  of  dust 
concentration;  the  increasing  heat  content  of  the 
lowest  levels  appears  to  b<  quickly  transferred  to 
the  mid-levels.  This  agrees  with  previous  modeling 


by  various  researchers  who  determined  that  the 
Martian  dust  storms  carry  the  seeds  of  their  own 
destruction;  increasing  the  static  stability  of  the 
lower  atmosphere  by  warming  middle  layers  inhibits 
the  dynamic  mechanisms  required  for  dust  lifting, 
and  the  storm  is  gradually  extinguished  from  below. 

Because  of  the  tendency  of  high  dust  loads  to 
decrease  the  amount  of  planetary  flux  reaching  the 
uppermost  layers,  the  maximum  dust  temperatures  and 
dust-gas  temperature  differentials  which  occur  at 
the  highest  levels  are  lower  for  the  thicker  dust 
loads.  Maximum  dust  temperatures  for  all  three 
simulated  dust  loadings  exceeded  280  K  with  an 
assumed  dust  emissivity  of  0.80,  while  maximum  dust- 
atmosphere  temperature  differences  exceeded  105  K 
at  the  42  km  level  in  each  case  during  daylight 
hours.  At  night,  the  minimum  temperature 
differentials  at  the  42  km  level  were  over  41  K  in 
each  case,  and  perceptible  dust -atmosphere 
temperature  differences  occured  at  lower  levels 
than  in  the  daytime. 

The  plotted  atmospheric  temperature  profiles 
indicate  that  at  least  five  distinct  atmospheric 
layers  exist  whenever  dust  is  present.  These 
include  a  surface  boundary  layer  of  about  1  km 
thickness  exhibiting  conduction/convection  dynamics 
and  experiencing  diurnal  temperature  shifts  of 


about  50  K;  a  well-mixea  region,  8  to  9  km  tnicv., 
of  adiabatic  lapse  rates,  just  above  the  boundary 
layer,  where  turbulence,  eddy  mixing,  convection, 
and  dust  influences  dominate;  a  30  km  thick  layer 
topping  the  adiabatic  region,  characterized  by 
radiative  heating/cooling  processes  with  diurnal 
temperature  changes  of  less  than  5  K,  and  whose 
static  stability  increases  with  increasing  dust 
load;  a  shallow  (2-3  km  thick)  layer  of  near- 
adiabatic  to  adiabatic  lapse  rates  overlying  the 
topmost  dust  level,  indicative  of  turbulent  mixing; 
and  finally,  a  region  which  returns  to  radiative 
heating/cool ing  processes.  In  the  dust-free  state, 
only  three  regimes  are  evident,  namely  the  surface 
boundary  layer,  the  well-mixed  adiabatic  layer 
above  it,  and  the  radiatively  active  layer 
overlying  them  both  and  extending  to  space.  It  is 
apparent  that  dust  loading  can  cause  significant 
changes  in  atmospheric  stratification  and  static 
stability;  detailed  discussion  of  the  attendant 
implications  for  the  pertinent  meteorological 
dynamics  (e.g. ,  storm  genesis  and  decay)  is  left 
for  future  study.  It  is  important  to  note,  however, 
that  on  Wars,  the  characteristic  relaxation  times 
for  thermal  perturbations  caused  by  radiation,  with 
spatial  scales  on  the  order  of  several  hundred 
meters,  are  comparable  to  the  relaxation  times  for 


thermal  perturbations  on  similar  spatial  scales  due 
to  mild  turbulent  convection  (Dannevick  arid 
Pallmann,  1974;  Belton  and  Goody,  1967);  this  is 
very  different  from  the  terrestrial  situation,  and 
suggests  that  the  character  of  dynamical  processes 
on  the  two  planets  can  be  quite  dissimilar. 
Consideration  of  the  tenuity  of  the  Martian 
atmosphere  also  supports  this  conclusion. 

Model  output  indicates  that  reflected  solar 
flux  at  any  level  is  about  50$  to  75$  of  the 
upwelling  ground  flux,  while  direct  incoming  solar 
flux  is  typically  more  than  1.5  to  2.5  times  as 
great  as  the  ground  flux,  depending  on  altitude. 
Atmospheric  layer  fluxes  ^re  about  20$  of  ground 
fluxes  at  a  given  level.  This  implies  that  solar 
fluxes  are  nominally  at  least  an  order  of  magnitude 
larger  than  atmospheric  fluxes.  Solar  flux  is 
greatest  at  the  highest  levels,  of  course,  while 
the  opposite  is  true  of  the  planetary  fluxes. 

Fluxes  due  to  particles  alone  depend  on  altitude 
and  dust  loading,  but  range  from  about  one  third  to 
twice  as  great  as  atmospheric  layer  fluxes.  It  is 
apparent  that  particle  flux  can  be  significant  when 
compared  with  the  atmospheric  flux,  but  since 
ground  emissions  dominate  the  amount  of  energy 
radiated  to  space  (output)  and  solar  energy  is  the 
primary  forcing  function  (input;  for  the  system. 


only  unusually  heavy  oust  concentrations  are 
capable  of  perturbing  the  overall  planetary  thermal 
balance. 

4.2  Validation  of  Modal  Output 

Volume  84  of  the  Journal  of  Geophysical 
Research  (JGR),  dated  June  10,  1979  and  published 
as  the  "Mars  Volatiles  Issue",  presented 
contributions  from  many  of  the  principal 
investigators  of  the  Viking  mission,  including 
several  articles  directly  pertaining  to  the  thermal 
infrared  properties  of  t:.e  Martian  atmosphere.  Dust 
aerosols  were  also  discusses  in  some  detail,  and 
derived  particle  parameters  were  incorporated  into 
the  model.  The  analyses  offered  in  the  Mars 
Volatiles  Issue  of  J'Jh  w.ll  be  usee  here  as  one 
benchmark  in  determining  trie  validity  of  the 
MARSDUOT  monel ,  since  they  represent  the  best 
understanding  of  cm  y.  .  ,n  atmosphere  yet 
available. 

The  previous .y  i  sc  rosed  attuner..ent  of  the 
monel  to  cyclical  balance  is  also  i  kind  of 
validation.  Obviously,  if  the  model  were  lr capable 
of  such  balance,  it  would  nut  even  be  internally 
consistent,  much  lens  ,  val.d  simulation  of 
reality.  In  addition,  a  necessary  validation 


procedure  Lo  to  compare  the  fluxes  computed  for  the 

dust  free  case  with  fluxes  computed  for  the  various 

dust  loads,  to  determine  if  the  "dusty"  atmospheric 

fluxes  are  approximately  equal  to  dust  free 

atmospheric  fluxes  plus  dust  particle  fluxes  alone. 

This  approach  is  more  susceptible  to  error  because 

of  the  lack  of  accurate  estimates  of  dust 

concentration,  which  are  needed  to  calculate  dust 

particle  fluxes,  but  it  is  hoped  that  even  a  rough 

order  of  magnitude  agreement  will  provide  increased 

confidence  in  the  model. 

In  this  regard,  upward  and  downward  fluxes 

computed  for  the  atmospheric  source  layers  by  the 

model  were  compared  for  the  dust  free  and  dust 

laden  cases.  Then  fluxes  due  to  dust  particles 

alone  were  computed  manually,  using  the  vertically 

integrated  column  number  density  derived  by  Conrath 

et  al.  (1973)  from  Mariner  9  data.  This  number 

density,  corresponding  to  the  MARSDUST  amplitude  of 
-7  -1 

1.2  x  1C  cm  ,  was  determined  to  be  about  1.2  x 

7 

10  for  a  radius  of  1  Am.  Assuming  an  atmospheric 

column  of  about  40  km  maximum  height,  the  Conrath 

et  al.  value  yields  a  mean  number  density  of  3/cm. 

Tnis  is  net  in  agreement  with  calculations  made 

later  by  Pollack  et  al.  (1977),  who  derived  a  value 
■? 

closer  to  30/cm  ;  however,  Pollack  et  al .  assumed  a 
radius  almost  an  order  of  magnitude  smaller,  which 
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probably  accounts  for  the  discrepancy.  1'  is 
important  to  note  that  these  uerived  number 
densities  apply  only  to  the  particle  radii  quotes, 
and  are  not  meant  to  describe  the  entire  size 
distribution.  The  Conratn  et  al.  figure,  being  the 
result  of  a  vertical  integration,  must  be  thought 
of  as  an  average,  implying  smaller  values  at  high 
altitudes  and  greater1  values  near  the  surface.  For 

3 

the  sake  of  this  analysis,  a  value  of  5/cm  was 

assumed  near  the  surface  for  a  DUSTAMP  (extinction 

coefficient;  of  1.2  x  10“"^  cm-'1';  similarly,  DUSTAMP 

-7  -1 

values  of  p.O  and  c.O  x  10  cm  correspond  to 

3 

near-surface  number  densities  of  about  12/cm  and 
25/cm-5,  respectively. 

In  a  rectangular  column  with  a  base  area  of 
one  square  centimeter  and  height  one  kilometer  (a 
typical  source  layer  thickness  in  the  model),  the 
total  number  of  particles  would  be  about  500  000, 

1  200  000,  and  2  500  000  for  the  three  DUSTAMPs. 
Since  tne  particles  arc  vastly  smaller  than  the 
volume  they  occupy,  ever,  these  numbers  are 
considered  to  induce  a  negligible  amount  of 
overlap  in  cross  section,  so  that  the  particles  are 
assumsc  not  to  obscure  ->r  "shade w"  each  other. 

Thus,  tne  principle  of  linear  superposition  can  be 
applied  to  compute  t:  «  net  upward  and  downward 
fluxes  from  such  an  ensemble  of  particles.  Using 


the  Stefari-Boltzmann  Law  fur  gray  bodies  and  an 
assumed  emissivity  of  0.80,  flux  densities  for 
individual  particles  at  a  given  temperature  can  be 
computed.  Knowing  their  total  surface  area  permits 
direct  calculation  of  the  total  flux  due  to  the 
particle  ensemble. 

Using  a  median  temperature  value  of  230  K  for 

a  1  km  thick  source  layer  centered  at  1.5  km 

altitude  at  local  noon  (12  LM3T),  the  sum  of  the 

total  upward  and  total  downward  particle  fluxes, 

based  on  the  Stefan-Boltzman  Law,  is  about  8000 

3 

erg/sec  for  the  5/crrr  dust  load,  19  000  erg/ sec  for 
3 

the  12/cm  concentration,  and  38  000  erg/sec  for 
3 

the  25/cm  concentration.  For  this  same  layer  in 
the  dust  free  state,  the  sum  of  the  upward  and 
downward  fluxes  is  about  19  000  erg/sec.  For  the 
three  DUSTAMPs,  in  increasing  order,  the  amount  of 
flux  in  excess  of  this  baseline  figure  for  the 
given  layer  is  10  000,  24  000,  and  43  000  erg/sec. 
These  values  differ  from  the  dust  fluxes  stated 
above  by  about  25$.  21$,  and  12$,  respectively.  In 
light  of  the  crude  assumptions  made,  this  agreement 
is  surprisingly  close  and  pennies  a  reasonable 
level  of  confidence  in  th '■  validity  of  the  model. 

Martin  et  al .  (1979)  gave  temperatures  for  the 
surface  and  dusty  atmosphere,  based  on  the 
inversion  of  Viking  IRTM  radiance  data  at  various 
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wavelengths,  which  are  in  agreement  with  MARSDU3T 
output.  They  conclude  that  surface  temperatures  are 
much  greater  than  atmospheric  temperatures 
typically  between  tne  hours  of  9  and  16  LMST,  a 
feature  captured  rather  well  by  MARSDUb'l .  According 
to  Martin  et  al.,  the  inversion  of  20  /im  radiance 
data  gives  an  estimate  of  minimum  surface 
temperatures  while  inverting  the  7  /m  data  yields 
maximum  surface  temperatures.  If  so,  their  dust- 
free  maximum  around  2l53  K  and  minimum  near  181  K  is 
in  good  agreement  with  the  model.  Their  dusty 
atmosphere  data  was  mainly  for  the  southern 
hemisphere,  where  summer  temperatures  are  warmer 
than  in  northern  summer,  but  even  so,  their  figures 
of  about  273  K  and  21u  K,  corresponding  roughly  to 
the  same  time  of  season  as  the  MARSDUST  simulation, 
are  very  close  to  model  outputs  for  the  heavier 
dust  loads.  Peterfreund  and  Kieffer  (19791 
determined  mean  dust  cloud  temperatures  ne  -r  237  K 
and  cloud  top  heights  between  8  and  17  km.  A 
derivation  of  dust  cloud  height  by  Martin  al. 
from  9  .m  radiance  data  implied  an  altitude  of  11 
km,  which  as  remarkably  close  tc  the  level  for 
minimum  dust  temperature  given  ;.y  the  model  f^r  the 
various  dust  loads.  For  model  output  belo*  this 
altitude,  adiabatic  laps*;  rales  ire  the  rub-,  while 
above  it,  slower  cooling  of  the  ntmospnene  w'.tn 


height  gives  increasing  stability  and  a  diminished 
capacity  to  raise  dust.  Indeed,  Thomas  and  Gierasch 
(1985)  report  that  local  "dust  devils"  rarely 
exceed  a  height  of  7  km-  Thus  the  level  of  minimum 
dust  temperature  can  indeed  be  considered  the  top 
of  a  dust  cloud  deck,  and  the  Viking  results 
support  this  interpretation.  The  amplitude  of 
diurnal  surface  temperature  changes  was  observed  by 
the  Viking  IRIK  to  be  within  a  range  of  about  50  to 
100  K,  while  typical  model  data  gives  amplitudes  on 
the  order  of  60  to  80  K. 

These  many  comparisons  are  offered  as  a 
validation  of  the  model;  in  no  respect  does  the 
MARSDUST  output  grossly  deviate  from  the  norms  set 
by  Viking  observations.  These  checks  of  internal 
consistency  and  similitude  of  output  to  actual 
measurements  and  known  phenomena  support  the 
contention  that  the  model  is  indeed  a  reasonable 


match  to  reality  on  Mars. 
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gas,  even  at  night.  For  the  Earth,  65  kilometers 
(corresponding  to  a  pressure  of  0.1  hPa  or  0.1  mb) 
is  the  threshold  altiude  above  which  this 
difference  becomes  larger  than  10  K,  increasing  to 
over  100  K  near  the  90  km  or  0.001  hPa  level 
(Fiocco  et  al.,  1975).  For  Mars,  the  current  model 
indicates  that  a  10  K  difference  can  be  achieved  at 
a  much  lower  altitude,  closer  to  15  km  (1  hPa), 
while  a  surplus  of  100  K  is  possible  well  below  50 
km  (indeed,  MARSDUST  restricts  dust  to  altitudes 
below  42  km,  or  the  0.1  hPa  level),  depending  on 
particle  concentration.  Because  of  the  higher 
density,  lower  infrared  emissivity,  and  more 
complicated  chemistry  of  the  terrestrial  atmosphere 
compared  to  that  of  Mars,  the  dependence  of  dust 
temperature  on  particle  radius  and  altitude  is 
stronger  and  more  complex  for  Earth  than  Mars.  Our 
model  neglects  this  weaker  dependence  of  dust 
temperature  on  particle  radius  for  Mars. 

The  atmospheric  temperature  profiles  derived 
from  MARSDUST  output  indicate  that  at  least  five 
distinct  atmospheric  layers  are  sustained  when  dust 
is  present:  a  surface  boundary  layer  of  about  1  km 
thickness  dominated  by  conduction/convection 
processes  and  experiencing  diurnal  temperature 
changes  on  the  order  of  50  K;  a  persistent  region 
of  adiabatic  lapse  rates,  8  to  9  km  thick,  just 


above  the  boundary  layer,  where  turbulent  mixing, 
convection,  and  dust  influences  are  maximized;  a 
layer  about  30  km  thick  above  the  adiabatic  region, 
characterized  by  radiative  heating/ coo ling 
processes  with  diurnal  temperature  changes  on  the 
order  of  5  K,  and  whose  static  stability  increases 
with  increasing  dust  load;  a  shallow  (2-3  km  thick) 
layer  of  near-adiabatic  to  adiabatic  lapse  rates 
overlying  the  topmost  dust  level;  and  finally,  a 
region  charaterized  by  a  return  to  radiative 
heating/cooling  processes.  Without  dust,  only  three 
domains  exist:  the  surface  boundary  layer,  the 
adiabatic  layer  above  it,  anu  the  region  of 
radiative  heating/cooling  overlying  them  both  and 
extending  to  space.  On  Mars  at  least,  it  appears 
that  dust  loading  can  cause  significant  changes  in 
atmospheric  stratification  and  static  stability. 

Atmospheric  temperatures  derived  by  inversion 
of  radiance  data  obtained  from  orbit  can  become 
suspect  when  aerosols  are  present.  Depending  on  the 
wavelength  channel  used  by  the  orbiting  radiometer, 
measured  upwelling  flux  can  be  substantially 
increased  by  net  particle  emissions,  which  may  be 
on  the  order  of  the  atmospheric  emission  itself  (or 
greater)  for  a  given  layer. 

Viking  Orbiter  Infrared  Thermal  Mapper  (IRTM) 
data  of  Mars  indicate  to  some  observers  that  sudden 


warmings  on  the  order  of  80  K  at  the  28  km  level 
near  the  poles  are  possible  and  are  linked  to 
atmospheric  dust  content  (Martin  and  Kieffer, 
1979).  One  may  wonder  if  this  is  an  exaggerated 


deposition  and  precipitation  of  CUi?  in  tin-  cold 
polar  atmosphere,  and  a  reduc'tion  in  the  dost 
cloud's  effective  emissivity  could  be  the  result,  of 
its  photon  energy  being  lost  to  slowing  the  phase 
change  process.  At  an  altitude  of  <?5  km,  oust 
concentration  (one  or  two  particles  of  1  /m  mean 
radius  per  err.')  and/or  atmospheric  density  (1.5  x 
10  ^  g/em^)  would  be  sufficiently  low  to  cast  doubt 
on  tne  conclusion  that  the  dust  is  directly 
responsible  for  any  significant  warming  of  the 
atmosphere  at  tnat  level,  iiowever,  an  unusually 
thiCK  but  necessarily  transient  cloud  of  limited 
vertical  extent  carried  poleward  by  Hadley  cell 
dynamics  (both  Earth  and  Mars  possess  atmospheric 
circulation  patterns  which  Lend  to  preferentially 
deliver  aerosol  particles  to  the  poles)  could  be 
the  source  of  significant  Ih  emission  before  it 
inevitably  succumbs  to  gravitational  settling  and 
polar  subsidence.  Inverted  temperatures  based  on 
radiance  data  under  such  circumstances  should  be 
viewed  with  a  higniy  critical  eye.  Our  model  output 
indicates  that  even  moderate  dust  loacs  may  greatly 
distort  the  nominal  flux  fields  one  would  expect  to 
observe . 

Kesesrcr.  utilizing  sa'cllite  radiometer  da ta 
should  tnerefore  be  sisenarveu  with  great  care, 
since  the  ev  i  a>'r, c  suggest:;  ft.. it  m  my  of  the  bi  t  -k 


body  assumptions  made  during  the  inverting  of 
radiance  data  can  be  invalidated  by  the  presence  of 
particulates,  especially  at  the  higher  altitudes, 
where  dust -gas  temperature  disparities  are  at  their 
extreme.  The  black  body  radiation  laws  are  the 
cornerstones  of  radiative  transfer  theory,  and  care 
must  be  taken  that  the  underlying  assumptions  of 
the  theory  are  not  violated  in  the  haste  to  arrive 
at  dramatic  conclusions. 

A  fully  time-dependent  planetary  radiative- 
convective  heat  transfer  model  is  beyond  present 
capabilities  because  the  complexities  of  orbital 
variations  and  intrinsic  planetary  structural 
changes  have  not  been  understood  or  described  for 
the  necessary  time  scales  (100  000  to  1  000  000 
years).  The  most  technically  advanced  models 
available,  such  as  the  Community  Climate  Models  of 
the  National  Center  for  Atmospheric  Research 
(NCAR),  attempt  to  address  these  issues  for  the 
terrestrial  case.  However,  restricting  model  output 
for  Mars  to  more  modest  human  time  scales  makes  the 
problem  more  tractable,  and  does  not  unduly  limit 
our  understanding. 

As  in  all  studies,  accurate  data  is  a  critical 
concern.  The  phenomenological  setting  for  Mars 
extends  us  to  the  limits  of  our  technology  with 
respect  to  data  acquisition  at  great  distances. 


This  significant  const  'air it  should  be  bonne  in  mind 
when  considering  experimental  procedures  and 
techniques  in  planetary  studies.  Much  desirable 
information  will  remain  unobtainable  for  a  long 
time  to  come.  Without  actual  dust,  soil,  or  rock 
samples  in  hand,  it  will  be  difficult  to  improve 
our  quantitative  understanding  of  the  spectroscopic 
properties  of  Martian  dust,  for  example.  Ihe 
insufficiency  of  the  observation  network  on  Earth 
is  still  hotly  debated;  it  is  not  hard  to  imagine 
how  our  understanding  of  any  other  planetary 
atmosphere  is  crippled  by  not  having  even  a  meager 
on-site  observation  capability  there. 

Despite  the  fact  that  disparate  planetary 
atmospheres  may  experience  similar  radiative 
effects  due  to  the  presence  of  particulates, 
differences  also  abound;  in  the  Martian  midlatitude 
atmosphere,  where  phase  changes  and  latent  heating 
are  negligible,  atmospheric  wanning  by  dust 
particles  can  act  to  replace  or  emulate  these 
processes.  Meridional  heat  transport,  storm 
genesis,  and  so  on,  do  not  seem  tc  require  the 
existence  of  volatile  components  within  tb-  Martian 
atmosphere,  unlike  the  terrestrial  case. 

Yet  the  extrapolation  of  terrestrial -based 
theory  to  other  planets  will  have  to  suffice  for 
the  irmieoiate  future;  this  approach  cm  provide  an 


adequate  foundation  for  further  research,  as  long 
as  tne  dissimilarities  between  differing 
atmospheres  are  kept  in  mind.  Extensive  modeling 
efforts,  such  as  MARSTEMP  (Pallmann  et  al.)  and 
LOWTRAM/FASCODE  (AFGL) ,  have  given  us  valuable 
insight  into  planetary  radiative  transfer.  The 
McClatchey  et  al.  and  Rothman  et  al.  compilations 
of  carbon  dioxide  spectroscopic  parameters  are 
evidence  of  what  Earth-bound  research  can  provide 
to  planetary  studies,  and  it  is  a  strong  step  in 
the  right  direction  for  Mars,  in  particular. 

But  for  the  Earth  itself,  we  lack  data  on 
basic  radiative  transfer  processes.  For  example, 
global  circulation  patterns,  cloud  effects,  aerosol 
loadings,  albedo  variations  (especially  for  ice), 
and  the  effects  of  radiatively  active  trace  gases 
are  all  imperfectly  understood  even  for  our  own 
planet — what  can  we  hope  to  say  about  these 
processes  on  stranger,  less  familiar,  and  far  more 


distant  worlds? 


APPENDIX 


VERTICAL  PROFILES  OF  ATMOSPHERIC  TEMPERATURE, 
DUST  AMPLITUDE,  DUST  TEMPERATURE,  AND 
DUST  TEMPERATURE  MINUS  ATMOSPHERIC  TEMPERATURE 

PLOTTED  FROM  SIMULATION  OUTPUT 
COMPUTED  FOR  VARIOUS  DUST  LOADINGS 


BY  THE  MARSDUST  PROGRAM 


TABLE  1.  PROGRAM  PARAMETERS; 


NO  DUST  LOAD 


Initial  Time:  16.000  LMST;  Local  Martian  Solar  Time 

Number  of  Timesteps:  1710 

Timestep  Interval:  10  minutes  (Martian) 

[4  minutes  (Martian)  from  LMST  13.500  to  16.500] 

Surface  Albedo:  0.284 
Surface  Emissivity:  0.900 
Soil  Density:  1.500  g/cur> 

Soil  Heat  Capacity:  1.01  x  10^  J/(m^  K) 

Soil  Conductivity:  0.230  J/(s  m  K) 

C02  Conductivity:  0.0133  J/(s  m  K) 

Latitude:  22.50°  N 
Solar  Declination:  24.60° 

Ratio  of  Mean  to  Actual  Solar  Distance:  0.9222 
(This  data  for  Viking  1  Lander;  20  Jul  1976) 


Surface  Pressure:  7.6  hPa  (mb) 

Dust  Extinction  Coefficient:  0.0 

One  Martian  Solar  Day  =  24  Martian  Hours  =  1440 
Martian  Minutes;  One  Martian  Minute  =  61.6625  s 
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MARTIAN  DUSTY  A'l'MOSF'Ht-RE  HAD-COND-OONV  HI' at 


TRANSFER  SIMULATION:  MARSDU3T 


Figures  A.l  -  A. 7 


Vertical  Profiles  of 
Atmospheric  Temperature 
With  No  Dust  Loading 
DUSTAMP  =0.0 


Four  hour  intervals  from  LMST  0.000  to  LMST  24.000 
Horizontal  axis  tempt; nature  range:  160  K  to  260  K; 
Vertical  axis  altituue  range:  C  m  to  50  000  m. 
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Atmospheric  Temperature  Profile;  No  Dust  Lose ,  LMST  24. 00 


TABLE  2.  PROGRAM  PARAMHL'ERG; 

DUSTAMP  =  1.2  x  10"7  cm-1 

Initial  Time:  16.000  LMST;  Local  Martian  Solar  Time 

Number  of  Timesteps:  676 

Timestep  Interval:  lu  minutes  (Martian) 

[4  minutes  (Martian)  from  LMST  13.500  to  lb. 500] 

Surface  Albedo:  0.284 
Surface  Elnissivity:  0.900 
Soil  Density:  1.500  g/cm^ 

Soil  Heat  Capacity:  1.01  x  10^  J/(m^  K) 

Soil  Conductivity:  0.2^0  J/(s  m  K) 

C02  Conductivity:  0.0133  J/(s  m  K) 

Latitude:  22.50°  N 
Solar  Declination:  24.60° 

Ratio  of  Mean  to  Actual  Solar  Distance:  0.9222 
(This  data  for  Viking  1  Lander:  20  Jul  1976 ) 

Surface  Pressure:  7.6  nPa  (mb) 

Dust  Extinction  Coefficient:  1.2  x  10-^  cm-'1' 

One  Martian  Solar  Day  =  24  Martian  Hours  =  1440 
Martian  Minutes;  One  Martian  Minute  =  61.6626  s 


MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HEAT 
TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 8  -  A. 18 

Vertical  Profiles  of 
Atmospheric  Temperature 
With  Light  Dust  Loading 
DUSTAMP  =  1.2  x  10-7  on-' 


Four  hour  intervals  from  LMST  20.000  to  LMST  4.000; 


two  hour  intervals  otherwise. 

Horizontal  axis  temperature  range:  160  K  to  260  K; 
Vertical  axis  altitude  range:  0  m  to  00  000  m. 
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Atmospheric  Temperature  Profile:  Light  Dust  lobo,  lMST  04.00 


mospheric  Temperature  Profile;  Light  Dust  Load,  lMST  15. 


i09pr>er>c  Temperst-jr^  Profile-  Light  Dust  Loaa  . 


MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HEAT 


TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 19  -  A. 29 


Vertical  Profiles  of 
Dust  Temperature 
With  Light  Dust  Loading 
DUSTAMP  =  1.2  x  lCf7  cnT 


Four  hour  intervals  from  LMST  20.000  to  LMST  4.000; 

two  hour  intervals  otherwise. 

Horizontal  axis  temperature  range:  190  K  to  290  K; 
Vertical  axis  altitude  range:  0  m  to  50  000  m. 


.19  Dust  Temperature  Profile;  L'otit  Dust  loafl ,  lMST  00.00 
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MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HEAT 


TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 30  -  A.AO 


Vertical  Profiles  of 

Dust  Temperature  Minus  Atmospheric  Temperature 
With  Light  Dust  Loading 
DUSTAMP  =  1.2  x  10-7  cnf1 


Four  hour  intervals  from  LMST  20.000  to  LMST  A. 000 
two  hour  intervals  otherwise. 

Horizontal  axis  temperature  difference  range: 

0  K  to  150  K; 

Vertical  axis  altitude  range:  0  m  to  50  000  m. 
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TABLE  3.  PROGRAM  PARAMETERS; 

DUSTAMP  =  3.0  x  10~7  cm-1 

Initial  Time:  16.000  LMST;  Local  Martian  Solar  Time 

Number  of  Timesteps:  576 

Timestep  Interval:  10  minutes  (Martian) 

[4  minutes  (Martian)  from  LMST  13.500  to  16.500] 

Surface  Albedo:  0.284 
Surface  Elmissivity:  0.900 
Soil  Density:  1.500  g/cm^ 

Soil  Heat  Capacity:  1.01  x  10^  J/(m^  K) 

Soil  Conductivity:  0.230  J/(s  m  K) 

C02  Conductivity:  0.0133  J/(s  m  K) 

Latitude:  22.50°  N 
Solar  Declination:  24.60° 

Ratio  of  Mean  to  Actual  Solar  Distance:  0.9222 
(This  data  for  Viking  1  Lander;  20  Jul  1976) 

Surface  Pressure:  7.6  hPa  (mb) 

-7  -1 

Dust  Extinction  Coefficient:  3.0  x  10  cm 

One  Martian  Solar  Day  =  24  Martian  Hours  =  1440 
Martian  Minutes;  One  Martian  Minute  =  61.6625  s 
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MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HKA' 


TRANSFER  SIMULATION:  MARSDUST 


Figures  A.4l  -  A. 51 


Vertical  Profiles  of 
Dust  Temperature 
With  Medium  Dust  Loading 
DUSTAMP  =  3*0  x  lu'7  cm' 


Four  hour  intervals  from  LM31  PO.OjO  to  LYST  4. 

two  hour  intervals  otherwise. 
Horizontal  axis  temperature  range:  lau  K  r.  o 
Vertical  axis  altitude  range:  0  rr.  to  LU  uf u  rr. 
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MARTIAN  DUSTY  ATMOSPHLRE  HAD -COND -CONY  HLAT 


TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 52  -  A. 62 


Vertical  Profiles  of 

Dust  Temperature  Minus  Atmospheric  Temperature 
With  Medium  Dust  Loading 
DUSTAMP  =  3.0  x  1G"7  cm"1 


Fear  hour  intervals  from  LMST  20.000  to  LMST  A. 000 
two  hour  intervals  otherwise. 

Horizontal  axis  temperature  difference  range: 

0  K  to  150  K; 

Vertical  axis  altitude  range:  0  m  to  50  000  m. 
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TABLE  4.  PROGRAM  PARAMETERS; 

DUSTAMP  =  6.0  x  10-7  cm-1 

Initial  Time:  16.000  LMST;  Local  Martian  Solar  T 

Number  of  Timesteps:  676 

Timestep  Interval:  10  minutes  (Martian) 

[4  minutes  (Martian)  from  LMST  13.500  to  16.500] 

Surface  Albedo:  0.284 

Surface  Emissivity:  0.900 

Soil  Density:  1.500  g/crrr1 

Soil  Heat  Capacity:  1.01  x  10^  J/(m^  K) 

Soil  Conductivity:  0.230  J/(s  m  K) 

C02  Conductivity:  0.0133  J/(s  m  K) 

Latitude:  22.50°  N 
Solar  Declination:  24.60° 

Ratio  of  Mean  to  Actual  Solar  Distance:  0.9222 
(This  data  for  Viking  1  Lander;  20  Jul  1976) 

Surface  Pressure:  7.6  hPa  (mb) 

-7  -1 

Dust  Extinction  Coefficient:  b.O  x  iO  cm 

One  Martian  Solar  Day  =  24  Martian  Hour's  =  1440 
Marti  in  Minutes;  One  Martian  Minute  =  01.6625  s 
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TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 63  -  A. 73 


Vertical  Profiles  of 
Atmospheric  Temperature 
With  Heavy  Dust  Loading 
DUSTAMP  =  6.0  x  10-7  cm- 


Four  hour  intervals  from  LMST  20.000  to  LMST  4.000 
two  hour  intervals  otherwise. 

Horizontal  axis  temperature  range:  160  K  to  260  K; 
Vertical  axis  altituoe  range:  0  m  to  50  000  m. 


63  A  t'Tios  pne  r  i  c  Temperature  Profile;  Heavy  Dust  Loao  ,  LMS .  00.00 
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emperature  Profile;  Heavy  Oust 


MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HEAT 
TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 74  -  A. 84 

Vertical  Profiles  of 
Dust  Temperature 
With  Heavy  Dust  Loading 
DUSTAMP  =  6.0  x  lCf7  cm" 


Four  hour  intervals  from  LMST  20.000  to  LMST  4.000; 
two  hour  intervals  otherwise. 


Horizontal  axis  temperature  range:  190  K  to  290  K; 
Vertical  axis  altitude  range:  0  m  to  50  000  m. 
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82  Oust  Temperature  Profile;  Heavy  Oust  Load.  LMST  18.00 
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MARTIAN  DUSTY  ATMOSPHERE  RAD-COND-CONV  HEAT 
TRANSFER  SIMULATION:  MARSDUST 


Figures  A. 85  -  A. 95 

Vertical  Profiles  of 

Dust  Temperature  Minus  Atmospheric  Temperature 
With  Heavy  Dust  Loading 
DUSTAMP  =  6.0  x  10-7  cm-1 


Four  hour  intervals  from  LMST  20.000  to  LMST  <4.000; 
two  hour  intervals  otherwise. 

Horizontal  axis  temperature  difference  range: 

0  K  to  150  K; 

Vertical  axis  altitude  range:  0  m  to  50  000  m. 


emoe rature  Prof i  'p :  He  a  v  .  Dust  Loac 


emoer a  tor e  Profile; 
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